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THE  EXISTENCE  OF  CONSERVATION  LAWS 


By 

Howard  Osborn 

Introduction 


> first  order  quasi-linear  homogeneous  partial  differential 
equation  in  the  independent  variables  x1,  i = l,...,n  and  the 
dependent  variables  u^,  j ~ l,...,p  is  a conservation  law  if  it 
is  of  the  form 


(C.l) 


0 


for  some  0^  which  are  f -unctions  of  x = (x^|  and  u = (u^j. 

If  a scalar  density  f = f (x)  is  fixed  on  the  space  of  x,  that  is, 
if  a change  of  variables  x1  = x'(x)  involves  p (x1)  = p (x) 


where  — ^ 

*x] 

(0.1)  as 


6 x 


*4*  i. 

is  the  Jacobian,  then  one  may  set  p = pT  and  write 


(0.2) 


l (-'iiC-l-i  ^0-1  ) = divT  = 0, 

o N v I r>  * * 


bx"  P 

which  states  that  the  vector  A with  components  "h1  is  solenoidal 
with  respect  to  the  density  P . (pThe  notation  div  is  useful  for 
changes  in  density;  for  example,  if  A is  a scalar,  then 


div  A = A div  A A. \ 
Ap  a p ’ 


I f a metric  ds  = 


^ gi^(x)dx1dx'- 


ij 


is  attached  to  the  space,  the  vector  A is  solenoidal  with  respect  to 
the  given  metric  when  p = J~ 


det  g.  . 
ij 


is  taken  as  the  scalar  density. 


In  order  to  investigate  whether  a given  equation  has  the  form 
(0.2),  or  in  order  to  investigate  whether  some  linear  combination  of 
several  equations  has  this  form,  which  is  the  purpose  of  this  paper, 


0i 


X, 


2- 


it  is  clearly  not  necessary  to  know  the  density  p . For  if  one  knows 
only  that  the  equation  is  in  the  form  of  a conservation  law  (0.1), 
then  for  any  density  ^ the  vector  'H'  ® ^ 0 is  solenoidal  with  respect 
to  that  density.  Henceforth  a density  will  be  introduced  only  for 
purposes  of  illustration. 

Conservation  laws  arise  in  several  ways  in  the  mechanics  of 
continuous  media.  If  a fluid  of  density  has  the  velocity  vector 
u - u(x)  for  x in  some  closed  bounded  region  R,  then  Gauss'  theorem 
asserts  that 


<0-’>  jj"r.Pd<r  * JJTd|V  " ?dX’ 


<b  R 

where  u is  the  component  of  u normal  to  the  boundary  'b  R of  R.  The 

left-hand  side  of  (0.3)  is  the  variation  of  the  total  mass  over  R. 

This  is  physically  clear  since  H gives  the  normal  rate  of  flow  throughout 

the  entire  boundary  bR.  If  the  flow  is  stationary,  then  clearly  the  total 

mass  over  any  R is  conserved,  so  that  the  integrand  on  the  right-hand  side 

of  (0.3)  vanishes,  showing  that  u is  solenoidal  with  respect  to  p . 

Variational  principles  furnish  another  source  of  conservation  laws. 

Let  f be  a function  of  u^.  u^,  and  x^" , where  u^  *■  i ” l,...,n, 

1 r 

j - l,...,p,  and  suppose  that  the  integral  J dx,  for  which  a stationary- 
value  is  sought,  is  invariant  under  a group  with  p parameters.  Then, 
according  to  a well-known  result  of  Emmy  Noether,  one  can  find  ^ linear 
combinations  of  the  p resulting  variational  equations  which  have  the  form 
of  conservation  laws  [ U ] . As  an  example,  if  f does  not  depend  explicitly 
on  the  independent  variables  x1,  then  these  variables  themselves  may  be 


.V 


-3- 


taken  as  uarameters.  Indeed,  by  appropriate  linear  combination  of 
the  variational  equations 


(O.U)  f . - 5*  — ^ — r f . = 0 

$1  bx1  UV 


one  finds 


(o.5)  -3  ♦ ...♦ 

>x 


«xx 


= C. 


(k-1, . . . ,n) 


where 


(0.6) 


01 

k 


r 

£ 

j-i 


Unfortunately,  the  system  (0.5)  is  not  necessarily  linearly  equivalent 
to  (O.U) , for  example  when  n < p. 

Occasionally  well-known  systems  of  equations  which  are  not 
customarily  derived  from  variational  principles  can  be  written  in 
the  fora  of  a system  of  conservation  laws.  For  example,  the  equations 


of  the  one-dimensional,  non-isentropic , non-steady  flow  of  an  adiabatic 
fluid,  where  p is  the  density,  u the  velocity,  p the  pressure,  and  S 
the  adiabatic  constant,  can  be  written  in  the  form 


-h- 


(0.8) 


div 

(1, 

u) 

f 

( u , 

p 

div 

u 

r 

div 

(U2 

♦i 

e 


2 ^ 


0 


2* 

r=i 


u p 


0. 


If  an  arbitrary  system  of  first  order  quasi-linear  homogeneous 
partial  differential  equations  is  given,  one  might  ask  if  it  is 
equivalent  in  some  sense  to  a system  of  conservation  laws,  or  at  least 
how  many  conservation  laws  can  be  obtained  from  it  by  linear  combination, 
as  in  the  preceding  examples . This  question  will  be  stated  more 
specifically  in  Chapter  I and  answered  for  a small  class  of  equations 
in  Chapter  II;  it  is  answered  in  a very  general  way  in  Chapter  III. 

Chapter  I is  of  an  algebraic  nature,  and  reduces  the  problem  of 
the  existence  of  conservation  laws  to  that  of  solving  an  over-determined 
system  of  linear  homogeneous  partial  differential  equations.  An  existence 
theorem  is  given  for  the  over -determined  system  in  Chapter  II,  under  the 
restriction  that  certain  integrability  conditions  are  satisfied.  For 
the  more  general  problem  solved  in  Chapter  III  the  question  is  more 
conveniently  stated  in  terms  of  exterior  differential  forms  in  order 
that  the  existence  theorem  of  Cartan  and  Kahler  can  easily  be  applied. 

A brief  sketch  of  the  theory  of  systems  of  exterior  differential  forms 
has  been  included  for  convenience. 

The  principle  results  of  this  investigation  suggest  that  in  general 
one  should  expect  systems  of  partial  differential  equations  to  contain 
many  more  conservation  laws  than  appear  on  the  surface.  For  example, 
when  n **  2 and  p = 3 every  system  in  which  the  coefficients  depend 


'Z.. 
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exclusively  on  the  unknowns  possesses  at  least  a two-parameter  family 
of  conservation  laws.  Even  more  surprising,  for  a wide  variety  of 
cases  one  can  find  conservation  laws  which  depend  not  only  on  a small 
number  of  parameters,  but  also  on  a certain  number  of  functions  of  a 
single  variable;  in  particular,  when  the  integrability  conditions  of 
Chapter  II  are  satisfied,  one  may  assign  p arbitrary  functions  of  one 
variable  t,o  obtain  infinitely  many  equivalent  systems  of  conservation 
laws.  In  other  cases,  although  certain  functions  again  may  be  assigned 
arbitrarily,  equivalent  systems  of  conservation  laws  do  not  exist. 


All  definitions  will  bo  indicated  by  underlining  new  terms  as 
they  arise,  with  an  accompanying  explanation  in  the  text. 

The  Einstein  summation  convention  is  occasionally  used.  The 
dummy  indices  will  always  be  lower  case  Greek  letters  and  the  range 
of  summation  l,...,p.  There  is  no  summation  over  Latin  indices  except 
when  indicated  by  ^ . 


CHAPTER  I 
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DEFINITIONS  AND  STATEMENT  OF  THE  FROBLEU 


1.  Well-determined  systems 

This  paper  deals  with  formal  properties  of  a system  of  quasi- 
linear  homogeneous  first  order  partial  differential  equations 


(1.1) 


n,P 

<1 


ki 


1 , * • • ,m 


in  the  independent  variables  x = {x1]  and  unknowns  u **  {u^j , where 
c<ki  ^ oC^i  (x  u) 

A system  is  dete mined  if  m ” p and  equations  (1.1)  are  linearly 
independent  at  any  given  point  (xo,uq)  in  the  product  space  XX  U of 
the  spaces  of  x and  u.  Clearly  the  determination  of  a system  is 
independent  of  the  coordinates  of  U since  a change  of  coordinates 
u «*  u(v)  merely  turns  (1.1)  into  a system 


(1.2) 


where 


1 


h ■ 


ki  iu.j 

3 


0, 


since  the  matrix 


is  non-singular,  the  rank  of  the  np  X p matrix  ')  is  the  same 

as  that  of  (c*-^). 

Determined  systems  exist  for  which  she  corresponding  Cauchy 


% 


■ £ ft-  ii.* 


problem  would  be  rather  unnatural.  For  example,  if  n = p * 2 the 


system 

f iu1 

(1.3) 

/ iu1 

is  determined.  Its  solution,  ux  an  arbitrary  constant  and  u an 

1 2 

arbitrary  function  f(x  ,x  ) , is  uniquely  determined  by  the  unusual 
1 2 12 

Cauchy  data  u = constant  and  u •*  f(x  ,x  ).  Similarly  the  system 


(l.U) 


( — - o 
\ ix1 


^5x^ 


0 


can  be  completely  solved  merely  by  assigning  the  Cauchy  data 
u1  = f1(x2),  u2  ~ f2(x2),  where  f^  and  f2  are  arbitrary. 

let  * - M be  indeterminates , and  for  each  i let  (°^^")  be  a 
p X p matrix  with  complex  entries,  {i  = l,...,n).  Then  pathological 
examples  like  (1.3)  and  (l.U)  may  be  eliminated  by  requiring  that 


n 

(1.5)  det  ( 1 l ± 0 

V i = 1 

be  a non-zero  form  of  degree  p in  J . Any  system  (1.1)  satisfying 
(1.5)  is  called  well-determined . Clearly  this  property  is  independent 


of  the  coordinates  in  X as  well  as  in  U. 


8- 


The J ^ may  be  assigned  complex  values  which  vary  covariantly 

with  coordinate  changes  m X=  Fon  any  point,  (x  .u  ) . the  annihilator 

“ o ' o'  ' 

of  any  co-vector  (|  ) for  which  the  determinant  (1.5)  vanishes  is 

called  a characteristic  element  in  X at  (x^u^).  Thus  a system  is 

well-determined  if  and  only  if  not  every  element  is  characteristic. 

If  (1.5)  contains  a linear  factor,  the  coefficient  of}\ 

appearing  in  this  factor  may  be  taken  as  the  i^  direction  number 

of  a characteristic  direction  in  X at  (x^,uq).  Geometrically,  if 

there  exists  an  (n-2)  - parameter  family  of  characteristic  elements 

in  X at  (x  ,u  ) which  have  a direction  in  common,  then  the  direction 
o-o  ’ 

is  characteristic.  For  n = 2,  when  p linear  factors  always  exist, 
the  preceding  definitions  coincide,  and  furthermore,  for  any  characteristic 
direction  one  can  find  a linear  combination  of  the  equations  of  the 
system  of  partial  differential  equations  such  that  every  unknown  appearing 
in  it  is  differentiated  in  the  given  characteristic  direction.  For  n > 2 
such  linear  combinations  do  not  necessarily  exist. 


| 2 . Tangent  spaces,  equivalence  of  well-determined  systems 

Let  £u^  * (u  s x}  be  a class  of  differentiable  functions 
u of  X.  By  choosing  a set  [u^, . . . ,u  ^ of  functionally  independent 
u's  about  any  given  point  of  X this  class  becomes  a differentiable 
manifold,  also  denoted  by  {uj  , with  the  values  of  u-^,...,u  as  the 
local  coordinates.  VTe  note  that  th-  operator 


q 

■ir.i 


= 1 h 


bu . 

J 

i 


x 


d 

sir 


is  a linear  differential  mapping  of  the  class 


--9“ 


L 


of  all  differentiable  functions  on  the  manifold  |u^  into  the  complex 
number  field,  at  any  given  point  of  X.  Since  any  operation 

i 

on  is  also  a linear  differential  mapping,  these  mappings  form  a 

vector  space,  the  tangent  space  of  {u\  , spanned  by  - , . . . «— . Two 

b xx  o xn 

vectors  in  this  space  are  identical  if  and  only  if  they  furnish  tne 

same  mapping.  Clearly  the  space  is  of  dimension  q. 

An  alternate  definition  can  be  given  merely  by  replacing  "the  class 

of  all  differentiable  functions  on  the  manifold  [uj  n by  "the  class  [u]  " 

alone,  or  even  by  , . . . ,u  ^ and  considering  the  operators  -X-.  This 

± q bx1 

definition  has  the  advantage  that  it  defines  a tangent  space  to  the 

class  £u}  without  constructing  the  corresponding  manifold.  As  an  example, 

consider  the  class  /u}  of  all  solutions  u of  tH^  **  0.  Its  tangent 

b xx  1 n 

space  is  just  the  annihiiator  at  any  point  of  the  vector  (<*  , . . . ,<*  ) j; 
clearly  there  is  no  point  in  constructing  the  manifold  {u}  to  discover 
this . 


The  tangent  space  just  defined  is  distinct  from  the  tangent  space 
of  U , which  is  a p-aimerisional  space  spanned  by  the  mappings 

bup 

of  all  differentiable  functions  on  the  manifold  U,  independently  of  X. 


Let  ^(u1, . . . ,up)^  = |(u'1 , . . . ,up)  si}  be  a class  of  p- tuples  of 

functions  u^  of  X.  Suppose  there  are  q p-tuples  (u^, . . . ,uP) , . . . , (u^, . . 
such  that  the  rank  of  the  q * np  matrix 
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K 


Ox 


n 


sp 

<lx  1 


iu1  *uP 

q q 
n • * ‘ \ ' "T  * * * 

ox 


fcu1  *uP  , 

q q / 

• • • ' • • • / 

\ n \ n / 

* ” ox 


o x 


is  q;  then  the  q p-tuples  are  defined  to  be  independent.  Any  operator 

is  a linear  differential  mapping 


M-’--  v^) 

of  the  class  ^(ux, . . . ,uP)^  into  p-tuples  of  complex  numbers,  at  any  point 


of  X,  where  the  f*-'s  are  arbitrary.  As  before,  if  there  exist  exactly  q 
independent  p-tuples  these  maps  form  a q-dimensional  vector  space,  the 
tangent  space  of  f(u\..., 

Suppose  that  (1.1)  is  determined  and  continuous  in  X/  U in  a 
neighborhood  of  (xo,uq).  Then  the  p vectors  = (*^, . . . . . . jOC^11) 

k = l,...,p  are  independent  in  the  usual  sense  for  each  (x,u)  in  the 
neighborhood  and  span  a p-dimensional  subspace  V of  some  fixed  np-dimensional 
complex  vector  space. 

Assume  that  for  any  well-determined  system 


(1.6) 


$i*i 


lei  ^u^ 


ix1 


0, 


k = 


the  tangent  space  T of  the  class  of  solutions  ( (u'1', . . . ,uP)}  of  (1.6) 
is  isomorphic  to  V*,  the  dual  of  V.  Then  one  can  define  two  systems  to 
be  equivalent  if  the  tangent  spaces  of  their  classes  of  solutions  are 
isomorphic  to  each  other.  Clearly  this  implies  that  they  possess  the  same 


solutions 


_11. 


Without  using  the  above  assumption,  two  systems  (1.6)  will  be 
called  equivalent  if  the  corresponding  spaces  V>  are  isomorphic.  This 
means  simply  that  one  can  apply  a non-singular  linear  transformation, 
depending  on  X/  U,  to  obtain  the  equations  of  one  system  from  those  of 
the  other.  This  definition  is  less  satisfactory  than  the  one  just 
suggested  only  in  that  it  characterizes  V*  purely  formally. 


jT 

l 

t 


b 


§ 3.  Statement  of  the  problem 

If  a given  system  is  equivalent  to  a system  one  of  whose  equations 

happens  to  be  a conservation  law,  the  original  system  is  said  to  contain 

the  conservation  law.  Geometrically,  the  space  dual  to  the  vector 

consisting  of  the  coefficients  of  the  — ^ appearing  in  the  conservation 

1 2 

law  contains  V*.  As  an  example  suppose  n - 2 and  p - 1.  Then  if et  and 
are  differentiable  functions  of  u alone,  the  system 

. 1 


(1.7) 


^x 


♦ °c2  = o 

<)  x 


contains  infinitely  many  conservation  laws,  depending  on  a single 
arbitrary  function.  Namely. 

(1.8)  -X-  f oC1(v).f(v)dv  ♦ f °<.2(v)f(v)  dv  = 0, 

^ix1  ^)x 


where  f=  f(v)  is  the  arbitrary  function.  Note  that  the  tentative 
assumption  of  §2  always  holds  for  p * 1 so  that  the  equivalence  of  (1.7) 
and  (1.8)  could  be  taken  in  the  earlier  sense  if  desired. 

The  problem  is  to  find  how  many  conservation  laws  are  contained  in  a 


2? 

* 
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Bn 


V 

^ ,,  .A. 


given  system.  In  particular  it  is  of  interest  to  learn  when  a system 
is  equivalent  to  another  system  consisting  entirely  of  conservation 
laws,  (0.8)  for  example,  and  in  how  many  ways  such  a representation 
can  be  found.  One  surprising  result  is  tnat  even  if  those  representations 
which  can  be  obtained  from  one  another  by  linear  combination  with  constant 
coefficients  are  included  in  the  same  equivalence  class,  there  is  still 
a wide  variety  of  systems  which  can  be  represented  by  infinitely  many 
equivalence  classes. 

From  now  on  we  restrict  ourselves  to  two  independent  variables, 

n * 2,  except  as  noted.  The  number  of  unknowns  will  be  arbitrary, 

although  all  of  the  chief  difficulties  of  the  problem  occur  when  p = 3. 

ki 

We  further  assume  that  the  functions  oC  . are  independent  of  X.  Thus 

the  existence  of  conservation  laws  is  a purely  local  problem  in  U. 
ki 

The  functions  o<.  are  assumed  to  be  as  manv  times  differentiable  as 

j 

needed  in  some  neighborhood  of  u = u;  they  are  assumed  to  be  analytic 
in  Chapter  III.  Only  well-determined  systems  will  be  considered.  Finally, 
in  order  to  avoid  considering  several  exceptional  cases,  the  p 
characteristic  directions  in  X will  be  assumed  to  depend  continuously 
on  U and  to  be  distinct,  but  r.ct  necessarily  real,  at  (xo,uq),  except 
as  noted.  Thus  it  does  not  matter  whether  the  system  is  totally 
hyperbolic  or  not,  e.g. 

For  convenience,  specific  reference  to  uq  or  (xo,uq)  is  usually 
omitted..  All  the  results  of  this  paper  will  be  given  only  for  a 
neighborhood  of  this  point.  If  the  appropriate  conditions  are  satisfied 
in  a larger  portion  of  U or  of  Xx  U,  the  results  can  clearly  be  extended 
to  this  portion. 


" 3 


§ U.  The  corres  ponding  linear  equations 


Let  the  system 


(1.9)  + ^4=0,  k-l,...,p 

r ox  r 6x 

Ji  j[  f 

satisfy  the  assumptions  of  ^3  and  let  5 = (J  5 ^ ) 

annihilate  the  ith 

characteristic  direction  in  X,  X * p. 

The  matrix  (/  ^ c<  i 4/2^^  ) rank  P-1,  so  that  to  each 

characteristic  direction  in  X there  corresponds  a unique  characteristic 
direction  in  U,  denoted  by  'Vi  Wx'-’ll  ) which  is  defined  by 

(1.10)  ( 3 1 **■  p1  * I 2 **  p2  = °»  ^n0  summation  on  A ) . 

Since  the  characteristic  directions  in  X are  distinct  the  matrix 
jj*)  is  non-singular.  Setting 

(1.11)  |3“- 

equation  (1.10)  becomes 

(1.13)  £ 1 Pjf  * 0 (k  = 1, • ,p) 

(no  summation  on  & ), 


so  that  the  ratio  ft  ^ l ,3  jP  is  independent  of  k.  Furthermore,  for 
at  least  one  value  of  k not  both  and  vanish  s ince  (1.9) 

is  well-determined,  that  is, 


(l.llo  det  < 5 ! Aki  * 1 2 P™  )/°  • 

Let  ) be  any  member  of  V*  and  set 


(1.15)  u?  *=  ^ 0”  f 


so  that  in  the  new  coordinates  V # is  defined  by 
(1.16)  ^pkl  <^1  + ^ 2 = °»  k-l,...,p. 

Since  s is  well-defined  and  independent  of  k we  may  use 

the  elements 

P _1  —l 

<r  = (0,...,  o-1,...,0j0,...,  as  a basis  in  V*  with 

respect  to  the  new  coordinates,  where 


0.17)  <-=*6. 11 


„ ■*  0,  (no  summation  on,^  ). 


But  (1.13)  and  (1.17)  imply  that 

/j  | 

(1.18)  det  ( 1 2„  I - °,  <£  " 1^-jPj 

V *"i  °~2  ' 

which  is  merely  the  definition  of  V#  in  the  new  cccrd: 
Suppose  (1.19)  contains  a conservation  law 

, i^1  ,)/  . 0 

(ia9)  ^ 5?  W J7 


Setting 


•£*w*jW •y-i'tt-  - S il< 
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(1.20) 


7)  j a~ 

( jo  i 


this  may  be  vrritten  as 

(1.21)  o~1  Up  01  + <r2  u 0=0, 


and  we  may  speak  of  the  tangent  vectors  U^,  in  the  tangent  space  of  U. 
Then  the  elements  which  arc  dual  to  a space  consisting  of  the  single 
vector  (U-^0  ,...,11^0';  Uj0  , ...,Up0  ) must  contain  V#,  expressed 
in  the  new  coordinates.  It  will  suffice  to  look  at  the  basis  elements 
of  V*.  For  each./  we  have 

_ H i _ P 2 0 

(1.22)  ff—1U^0  2 (no  sumrna'tion  on  X , 

^ * l,...,p). 


Referring  to  (1.18)  this  becomes 


> i 1 Ar  p 

(1.23)  / 1 U^0X  ♦ 5 2 U 0‘ - - 0, 


(no  summation  on  X, 

X “ l,...,p) • 


It  should  be  noted  that  these  equations  are  entirely  independent  of  X. 
Except  for  the  trivial  constant  solutions,  there  corresponds  a conservation 
law  to  every  solution  of  (1,23) ? which  is  a linear  system  of  partial 
differential  equations  in  the  space  U. 
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5 , The  case  p = 2 

In  this  case  the  tangent  vectors  in  U are  of  the  form 

a.2M  • 

Let  (5  be  the  inverse  of  (?|k),  ^2^±  = <^±* 


Then  an  integrating  factor 


M. 

J 


always  exists  such  that 


Mj(  5 ^ du^"  * $ 2 ^u2)  1 total  differential,  say  dv^,  j *■  1,2, 


so  that 

(1.25) 


M P J - dv,“ 

j s k *7' 


(no  suinmation)  * 


hence 


(1.26) 


7 


k = M ^ U 


1 ^vJ 


and  (1.23)  becomes 


1 2 

(1.27)  $ i ^-r  + * 0,  j * 1,2,  after  division  by  M.. 

The  system  (1.27)  has  a number  of  solutions  about  uq  defined  by  two 
arbitrary  functions  of  a single  variable,  as  one  easily  shows  by  the 
Cauchy-Kowalewski  theorem,  if  the  } V are  assumed  to  be  analytic. 

This  existence  theorem  will  be  given  in  a particularly  nice  form  as  a 

special  case  of  a more  general  result  in  Chapter  III,  | §3,  5«  Thus 

for  p - 2 (1,9)  is  equivalent  to  many  systems  of  conservation  laws  [ .3  ] . 


-17- 


§ 6.  The  case  p > 2 

In  general  one  cannot  expect  to  find  integrating  factors  leading 
to  the  simple  form  (1.27)  of  equations  (1.23);  the  case  where  this 
can  be  done  is  treated  in  Chapter  II. 

In  any  case  the  system  (1.23),  or  any  system  equivalent  to  it,  is 
over-determined  for  p > 2,  that  is,  there  are  more  linearly  independent 
equations  than  unknowns.  Therefore,  in  order  to  find  results  analogous 
to  those  found  in  §3  and  £5  for  p * 1 and  p « 2 respectively,  one 
should  expect  that  certain  integrability  conditions  must  be  satisfied. 

As  noted  in  J U the  matrix  (0?^)  is  non-singular,  hence  the  p 
operators  are  differentiations  in  p distinct  directions  which 

span  the  tangent  space  in  any  point  of  U.  The  do  not  commute  with 

each  other  in  general.  However,  their  commutators 


(1.28)  [U^  Uffl] 


again  lie  in  the  tangent  space  and  so  may  be  expressed  in  the  form 


(1.29) 


U. 

.1 


for  certain  analytic  functions 

i.m  ira 

functions . One  easily  checks  that  for  any 
function  0 = 0 (u) 


(1.30)  u^CiyJ)  - um(y)  = [yjra]  0 


(u)  called  the  structure 
suitably  differentiable 


-18- 


Cleariy 


(1.31)  y J ♦ = o ; 

J!m  m!L 


furthermore  the  Jf J 's  must  satisfy  a relation  arising  from 
i.m 

Jacobi's  identity 

2 ^ufUJU  1 s 

t JL  m n 

Ann 

(1.32)  [[u^n)Dn]  • ttiy^v . nudity  =0, 


where  Z represents  cyclic  summation  as  indicated.  If  the  ]J  ^ 's  are 

constant,  the  fundamental  theorem  of  Lie  groups  asserts  that  there  exists 

an  analytic  group  for  which  the  are  the  infinitesimal  transformations. 

If  =0  except  when  i 'I  or  J - m one  can  find  integrating 

Im 

factors  M £ and  a change  of  variables,  v * V(u),  so  that 


(1.33)  U ^ s Uj  -M,  i=  1,...,P 


as  in  J5.  This  will  be  shown  in  Chapter  II,  which  deals  entirely  with 
this  special  case,  and  again  by  a simpler  means  in  Chapter  III. 


^ 7.  A preliminary  transformation 

The  elementary  methods  of  Chapter  II  are  better  suited  to  deal  with 

the  case  where  either  ( oC  ^ ) or  ( °t  . ) is  a non-singular  matrix. 

J J 

Suppose  that  ( oC  ) is  singular.  Since  (1.9)  is  well-determined  there 
J 
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exist  constants  - (/^.o  such  that  ( j*-  ^ °*-  ^ + **■  ^ ) 

1 2 

is  non-singular  in  a neighborhood  of  u “ uq  and  such  that  jU-  and  /*- 
are  distinct  directions  in  the  dual  of  the  tangent  space  of  X.  letting 


(1.3U)  x1  - x1  ♦ /*£  x2  , i - 1,2 


so  that 


(1.35)  -V 
bx1 


we  see  that  (1.9)  nay  be  written 

(1.36)  $ (^fczi  * * f ) 

J a x J lx 

1»1 


where 


(1.37)  - /*{<*■  f 


and  (<xk±)  is  non-singular.  Thus  (<*kl)  might  as  well  be  assumed 
J J 

non-singular  at  the  outset. 


(1.38)  H-  £ >?j  <r1 

*x  M 


(1.39)  <J-±  •/*■}  °'1  ^2. 
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w 

f. 


fc 


But  the  form 


_ 6 -|  0 
(1.U0)  o-  , U:  O1  ♦ <r~ 


oi  a conservation  law  is  invariant  with  respect  to  linear  coordinate 
clianges  in  X.  Hence  the  existence  conditions  ' (1.23)  become 


(1-1*1)  J-f  U^1  ♦ J g U#2  = 0 


jj  1 — Z 2 —H  r X 

where  / . » /<  ^ ♦ /*  -^2"  ■^1US  the  ratios  5^  s *2  an^ 

point  of  U may  be  changed  by  a projective  transformation.  In  particular, 

since  the  projective  group  is  simply  transitive  on  any  three  of  these  ratios, 

for  p - 3 one  may  adjust  the  characteristic  directions  in  X to  be  any 

three  arbitrary  distinct  directions  at  uq.  This  transformation  could 

have  been  derived  directly  from  (1.23)  simply  by  setting 


- I . 2 A 


(1.1*2)  ^ 01  ♦ /*-\  0' 


i nf2 


without  investigating  the  corresponding  transformation  in  X. 

Taking  ( oC  . ) to  be  non-singular  merely  corresponds  to 

n ‘ 2 


JL  = l,...,p. 


kl 

Assuming  that  ( . ) is  non-singular,  equations  (1.9)  can  be 

J 

written  in  the  unsymmetric  form 


(1.1*3)  ^ * 


hi  x ^ * k hi  = 0 
2 J 


j=i 


• • • y P* 


a 

-v 

vf 


1 
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In  this  case  the  derivation  of  § U is  easier  to  give.  For  suppose 
a conservation  law  exists,  so  that 

(1.19)  ^ (^i  W . 40!  w , .0, 

huj  ox  buj  bx 


Multiplying  (1.U3)  by 


and  adding, 

buK 


(l.UU)  £ 


buk 

*7  *7 


PJ-P  Wl 


< ^>0  o<  k 

.f  . buk  J 

J,k*JL 


hence 


/,  P<  b0X  * k buj  _ <'  b02  buj 

.7  . buk  - 6x^  .S-  buj  bx 

j,k*l  J~1 

But  the  p elements  ^=— *•,  j » l,...,p  determine  a basis  in  V#, 

bx^ 


so  that  (1.45)  implies 


(1.1.1)  i 


J = l».***p 


from  which  one  derives  (1.23)  merely  by  diagonalizing  ( °<  . ), 

J 

which  is  possible  because  the  characteristic  directions  in  X are  distinct. 


CHAPTER  II 


THE  SPECIAL  CASE 


§1,  The  tangent  vectors  IK 


A set  of  tangent  vectors  (operators)  U . , j * l,...,q  is 

complete  if  and  only  if  the  commutator  of  any  pair  of  them  lies  in 

the  tangent  space  the  set  spans.  (For  convenience  here  we  speak  of 

tangent  vectors  even  in  the  absence  of  a manifold  or  class  on  which  a 

tangent  vector  could  be  defined.)  It  is  well-known  that  for  any 

complete  set,  depending  continuously  on  r parameters,  r > q,  then 

there  exists  an  r - q parameter  family  of  integral  manifolds  in  the 

parameter  space  whose  tangent  spaces  are  spanned  by  the  given  vectors. 

Given  a set  of  vectors  U.,  j “ l,...,p,  suppose  there  exist 

I 

non-zero  factors  M.  and  a change  of  variables,  v « v(u),  such  that 
J 

these  vectors  are  of  the  form 


(2.1) 


U . £ M . 

J J bvj 


J = 1. 


»P- 


Then 


(2.2)  [U.U. ] = 5 1.  U.  - 5 U. 

v ] i ji  l ji  j 


where 


(2.3)  M.  t ±..  - U.M.  , 

i Ji  J i 


all  the  other  structure  functions  vanishing  identically.  This  is  clearly 

equivalent  to  the  assertion  that  any  subset  of  the  vectors  U.,  j = l,...,p, 

1 
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is  complete.  Conversely,  the  following  theorem  holds. 

THEOREM  At  If  any  subset  of  tne  vectors  U.,  j = l,...,p  is  complete, 

2 

then  there  exist  non-zero  factors  M.  and  a change  of  variables, 

2 ’ 

v = v(u),  satisfying  (2.1). 

Proof:  A sequence  of  changes  of  variables  will  be  constructed  in  such 

t h 

a way  that  at  the  n stage  the  form  (2.1)  is  displayed  for  all  j, 

< 

j = l,....n  n - p.  Assume  that  the  first  n - 1 stages  have  been 
completed.  Then,  after  multiplication  by  suitable  scalers,  the  first 

n vectors  may  be  assumed  to  be  of  the  form  n , 

p ou  6u 

"7  , where  - 1.  Since  by  assumption  the  commutator  ] 

du^  ^uj  duf 

is  a combination  of  — — and  we  have  "7  ^ “ A ^ ^or 

Ouj  n ouJ  J 

some  A , j / k.  But  ^ 11  - 1 implies  ^2-  - 0,  hence  A . * 0, 

J 6uj  J 

so  that  72^  - ^(uA?  un,...,uP)  for  j - l,...,n-l,  and 

'Y^  * ^ Au11, , . . ,uP)  for  j = n, . . . ,p. 
th 

The  n ‘ change  of  variables  will  be  assumed  to  have  the  form 


v2  . v')(u'-,3  un, ...  ,u^) 
vJ  - vJ(un,...,uP) 


(2.U) 


j - 1, . . . ,n-l 


j - n,...,p. 


Then  for  j ■ l,...,n-l 

(2.5)  U.  5 — e s , (no  summation  on  j) 

^ duj  duj  V?  duj  &vj 


so  that  the  first  n-1  vectors  retain  the  desired  form  provided 


I 
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f 0.  The  n^  vector  becomes 
<Kij 


(2.6)  U nrf  i—  s7i?  H.  -A_ 
n i/  ' iu?  ^v0"’ 


and  hence  will  have  the  form  U = ^AL.  — — if  and  only  if 

n ' *uP 

r(P  h:  A 


vanishes  for  j ^ n.  For  j * l,...,n-l  this  requirement  is 


(2.7) 


’f  i 


* S v k = °> 

^uj  k=]_  ^ 


j * 1, . . . >n— 1 , 


and  for  each  of  ohese  equations  one  can  find  a v^  for  which  / 0. 

&uj 

To  satisfy  the  requirement  for  j = n,...,p,  each  of  the  vJ,s,  j = n,...,p, 


may  be  chosen  as  an  independent  solution  of  the  equation 


(2.8) 


0, 


in  which  the  coefficients  depend  only  on  un,...,u^.  In  particular, 
vn  may  be  chosen  in  such  a way  that  f 0.  Clearly  the  Jacobian 

bun 

matrix  of  the  transformation  (2.1;)  chosen  in  this  way  is  non-singular. 

For  there  is  nothing  below  the  first  n - 1 terms  on  the  principal 
diagonal,  and  the  block  of  side  p - n ♦ 1 in  the  lower  right-hand 
corner  represents  the  ncn-singular  change  of  variables  chosen  last. 

This  completes  the  proof. 

Some  further  properties  of  vectors  satisfying  (2.1)  will  be  given 
for  convenience.  First,  for  any  differentiable  function  f,  (2.3)  implies 


J 


_oc__ 


(2.9) 


- U.f 
uj  (_m:)  = “m7 

i i 


U .M. 

f 

(M.)2 


* ~ (u.f  ♦ f jr1 . ) - i_  u. 

M.  ' J M.  o;  i 


where  U.  , is  the  operator  M,  U . (-rj—  ) s U.  ♦ S ^ . 

J5K  k j J ij 


If  f 


\ f 

i - — , j = then 

J d vj 

if.  'if.  , 

0 = — i i = U.f  . - ~ U.f. 

iv1  ivj  Ui  1 J Mj  ^ 1 


(2.10) 


i r-k-  fu.  . (U.f.)  - U.  . (M.f.)\ 
U7UT  (_  x,j  j y jii  i i'J 


In  particular,  since  (2.3)  can  be  written 
k 


<2’u)  ^ 


we  have 


(2.12) 


Ui;j*jk  Uj;iJikJ  i’"»k  " 1>***»P* 


Equation  (2.10)  can  be  given  an  interpretation  in  terms  of  the 

operators  U alone.  Given  p functions  g.,  a necessary  and  sufficient 
J 1 

condition  that  there  exist  a single  function  g such  that  g^  " U.g, 
j * 1, . . . , p is  that 


U.  .g.  = U.  .g., 
j;i  i i;j  j’ 


i j j * l,...,p. 


•i-7TvTr 


(2.13) 


'5,'*  p ^ ** >*pr,  _>  f :"-  *f*  »-r». 
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Another  proof  of  this  statement  follows  from  a form  of  Stokes’  theorem. 

let  Z represent  cyclic  summation,  Z F(ijk)  = F(ijk)  ♦ F(kij)  ♦ F(jki), 
ijk  ijk- 

i 0 — i 

and  suppose  0^  du  = u>  , then  for  p « 3 we  have 


(2.lU) 


f dg  = r (3  uj  ^ 

&R  on 


V du^du^ 


where  E * det  Formula  (2.lh)  can  be  proved  by  means  of  the  more 

customary  form  of  the  Stokes  theorem. 

The  factors  M.  can  be  calculated  directly  from  (2.3)  without 
J 

involving  the  sequence  of  changes  of  variable  used  in  the  proof  of 
THEOREM  A.  Once  these  factors  are  known  the  desired  change  of  variables 
can  be  made  directly. 

An  alternate  proof  of  THEOREM  A will  be  sketched,  based  on  the 

direct  computation  just  suggested.  By  performing  the  calculation  indicated 

in  Jacobi’s  identity,  (1.32).  substituting  tf^Up  for  [IKU^]  whenever 

v k 

possible,  one  finds  (2.12)  under  the  assumption  that  0 „ = 0 except 
for  k " j or  i = k.  For  each  i the  system  (2.3)  can  be  ccnverted 
into  a homogeneous  system  by  considering  as  a new  independent  variable; 

then  M4  can  be  found  in  terms  of  the  other  independent  variables,  if  the 
new  system  is  complete,  simply  by  considering  the  constant  solutions  to 
the  new  system.  It  turns  out  that  (2.12)  is  exactly  the  condition  that 
the  homogeneous  system  is  complete,  hence  the  factors  Kh  exist.  Thus 


i 
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the  commutators  of  the  vectors  - — U.  all  vanish,  from  which  it 

U . 

i 

follows  that  they  are  of  the  desired  form,  -^-r  . 

'brx 

A third  proof  of  THEOREM  A is  given  in  Chapter  III,  ^2. 

In  the  new  coordinates  the  axis  may  be  described  as  the 

solution  of 


(2.15) 


d / _ = dy^  = . dyP 

0 "■  1 "*  0 


passing  through  the  origin.  In  the  original  coordinates  this  is  the 
solution  of 


(2.16) 


du 


du* 


passing  through  uq.  The  curves  thus  obtained  for  j - l,...,p  are  at 

each  point  tangent  to  a characteristic  direction  in  U for  the  system 

(1.9)  and  are  called  the  characteristic  curves  through  uq.  As  THEOREM 

A has  shown,  it  is  the  peculiarity  of  the  special  case  being  considered 

in  this  chapter,  in  which  any  subset  of  the  tangent  vectors  U.  is 

J 

assumed  to  be  complete,  that  the  p families  K.  of  characteristic 

J 

curves  through  all  points  of  U may  be  used  to  define  a coordinate 
system.  This  means  that  if  .j-, , . ..,j„  is  any  permutation  of  l,...,p 

— i p 

then  those  members  of  K.  intersecting  given  members  of  K.  form 

d2  Jl 

the  same  family  K.  . of  two-dimensional  surfaces  when  j,  and  j » 

0^2  x 

are  interchanged;  those  members  of  K.  intersecting  given  members 

J3 


I 

i. 
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of  K.  . form  t.ie  same  family  of  three-dimensional  surfaces  when 
JjJ  2 

Jl  ^2  are  Permuted } etc.  These  hypersurfaces  are  called  the  webs 
spanned  by  j1  j2,  ^ etc. 


^ 2 . Conditions  for  a representation  by  a system  of  conservation  laws 

Assume  that  ar.  appropriate  change  of  variables  has  been  performed 

kl 

as  in  Chapter  I,  §7,  so  that  ( c<  ' ) is  non-singular.  Then 

JJ/O,  j =l,...,p  so  that  in  the  special  case,  in  which  any  subset 
of  the  tangent  vectors  U.  is  complete,  the  existence  conditions  (1.23) 
become 


(2.17)  = 0, 

bvj  bvj 


j * 3., . . . ,p 


after  division  by  M_.,  where  5 2 * Ji  / • Note  that  the  are 

distinct.  The  number  of  solutions  of  (2.17)  gives  the  number  of 

conservation  laws  contained  in  (1.9).  More  specifically,  the  question 

here  is  to  represent  (1.9)  by  a system  of  conservation  laws.  The 

necessary  and  sufficient  condition  is  that  there  exist  solutions 

( . . . ,^-j  ) cr  C2-1?)  which  span  V.  Hence, 

^v1  ivp  V ivp  2 

from  (2.17),  where  the  values  of  ^Cr,...,^-  determine  those  of 
wl  X(/l  ^ ^v-p 

5^-.....^--.  it  is  clear  that  a necessary  condition  is  that  no  linear 

^ Vp’  ,*2  ^ 

relation  exist  among  the  ^Tj***}!^* 

dy  0vp 

We  try  to  concoct  such  a linear  relation.  Let  i j k represent 
any  three  distinct  integers  among  i,...,p. 


b2#1  _ iV 

Since  ~ — r -tz 

0vK  ivj 


we  may  cross-differentiate  the  equations 
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of  (2.17)  to  find 


(2.ie)  Mi  y.  . }i  L i_  . «J_  fcl  . i * 1 i _ _ , 

ivk  ivj  <>vk&vj  ivj  )v:  £vk 


that  is,  setting  0 s 0 for  convenience, 

(2.19)  r * — — r - Xl  X_  ) 

*vj*r  X~?k  b v j dvk  ^vk  bvJ 


since  J ° / J k . It  should  be  noted  that  the  system  of  ( ^ ) second 
order  equations  obtained  by  writing  down  every  instance  of  (2.19)  is 
entirely  equivalent  to  (2.17).  Differentiating  (2.19)  with  respect 
to  v1  gives  rise  to  new  mixed  second  derivatives  on  the  right  hand 
side,  which  can  be  evaluated  exactly  as  in  (2.19)  in  terms  of  the  first 
derivatives  of  0 . Similarly  the  third  mixed  derivative  - 0 « . — c-  can 

be  expressed  in  terms  of  the  first  derivatives.  Equating  the  two  third 
derivatives  gives  a linear  relation  among  the  first  derivatives  of  0 


all  of  whose  coefficients  must  vanish  if  the  vectors 
1 2.  2 2 

(^t,...,^-  ; ) span  V.  The  coefficient  of  in  this 

bv1-  dr?  dvL  iv*5  bvK 

linear  relation  is  especially  easy  to  compute  if  one  notes  that  the  mixed 

2 *2 

second  derivative  - * - — - — r is  expressed  as  a linear 

bv1  ivj  bvj  bv1 

X0  A 0 

combination  of  only  the  first  derivatives  ~-r  and  . Equating  this 

o v1  bvj 

coefficient  to  zero  gives 

(2.20)  -X  (~— g ) - -X  (-X— ^ X-  ) - 0 

bv1  Jf J-  Jk  bvj  bvj  Ji-Jk  Av1 

as  the  necessary  condition. 
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Or.e  might  expect  that  the  coefficients  of 

lead  to  new  relations  which  cannot  be  obtained  merely  by  permuting  the 

indices  i j k in  (2.20).  They  do  not.  In  fact,  lecting  ©k . 

* 3 

represent  the  left-hand  side  of  (2.20),  the  entire  first  order  equation 
in  0 becomes 

0k 

(2.21)  Z 13  60  _ A 

ijk  ^vk 


- and  would 

'W1  bvj 


where  Z represents  cyclic  summation, 


Z F(ijk)  - F(ijk)  + F(jki)  ♦ F(kij). 


ijk 

Equation  (2,21)  can  be  proved  by  direct  calculation,  which  is  not  very 
instructive.  A simple  reason  that  one  set  of  indices  i j k leads  to 
only  one  linear  relation  among  the  first  derivatives  of  0 will  be  given 
in  Chapter  III,  §3. 

Equation  (2.21),  when  written  back  in  terms  of  the  original 

variables  and  tangent  vectors  tb , could  be  derived  directly  without 

making  a change  of  coordinates.  This  might  lead  one  to  expect  that  a 

similar  derivation  with  any  tangent  vectors  U.,  not  restricted  to  the 

3 

special  case,  would  give  a similar  result.  In  fact  an  analogue  to  (2.21) 
can  be  found,  but  unfortunately  it  is  no  longer  in  general  a first  order 


equation;  indeed  it  contains  terms 


J1-?! 


(Uj^0),  where  i j k are 


distinct.  Thus  the  only  case  in  which  the  resulting  equation  is  of 
first  order  is  the  special  case  considered  in  this  chapter.  The  more 
general  second  order  equation  can  of  course  be  used  to  derive  further 
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results  in  the  general  case,  however,  which  will  be  done  in  a different 
setting  in  Chapter  III. 

The  integrability  condition  (2.20)  for  the  special  case  will  be 
more  useful  expressed  in  terms  of  the  original  variables  and  tangent 
vectors,  since  the  coordinate  transformation,  v • v(u) , and  the 
integrating  factors  M might  be  inconvenient  to  find  in  any  given 
example . 


ar 


Setting  fk  = ^ ■ . 

3 Jj-Jk  *v.i 


, we  note  that  according 


to  (2.10)  relation  (2.20)  may  be  written  as 
(2.22)  U.  .(M.  f*)  - U.  .(M.  fk), 


that  is 


(2.23)  U,.,(— K 


u.r 

JL ) 

1J3  '!Mk 


u.lk 

uj;i  M ‘ykJ 


ix-y 


Thus  we  have  proved 


THEOREM  B;  Suppose  <*k  = k 


.(u)  is  a p X p matrix  of  twice 

J j 

differentiable  functions  of  u which  has  p distinct  characteristic 


roots 


Then  the  system 


(?.2h)  ♦ «<■  k =o  v = i ... 

• * ' ^ T er- l 7 


may  be  expressed  in  the  form 

. o — o- 

(2.25)  ^ ^ 1 j ^ =0,  k = l,...,p. 

&x  dxc 


Let  lb  represent  the  tangent  vector  ^ , 


P 1 


and  suppose  that 
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any  subset  of  the  vectors  U^,...,U  is  complete.  Then  a necessary 
condition  that  (2.2U)  is  equivalent  to  a system  of  conservation  laws 
is  that 

u.$k  u.yk 

(2.23)  u (-*-*)  - U (^-k) 

1,J  Jj-r  3 -s 

for  every  distinct  i j k among  l,...,p.  According  to  (2.13) 

1c 

it  rs  equivalent  to  demand  that  there  exist  p functions  f such 


(2. 214)  0,fk  - -i—., 
° })-} 


j,k  = l,...,p,  j / k. 


Clearly  any  system  with  constant  coefficients  satisfies  the 
conditions  of  the  preceding  theorem.  To  construct  a non-trivial  system 
satisfying  these  conditions  one  might  merely  write  down  a system  of 
conservation  laws,  attempting  to  choose  p pairs  of  functions  in  such 
a way  that  the  resulting  system  falls  into  the  special  case.  This  method 
is  much  more  difficult  than  attempting  to  guess  the  functions  J 1 directly. 
It  is  easy  to  check  that  the  following  example  satisfies  the  conditions: 


(2.25) 


, 23  ^u  _ _ 

=r  ♦ U U 5-  = 0 

^ x~  6x^ 


? 2 
bu  . 3 1 bu 

y "*■  U U " mK 


I ou  1. 

I T ♦ U U 

\ ^Xx 


12  bu 

u r? 

ox 


Another  simple  example  is  that  in  which  7^  = J^(u^),  i = l,...,p, 
whose  solution  one  can  easily  give  by  assigning  p arbitrary  functions 
exactly  as  one  was  assigned  in  (1.8). 


§ 3.  Proof  that  the  preceding  conditions  are  sufficient 

According  to  THEOREM  B there  exist  functions  f such  that  in 
an  appropriate  coordinate  system 


/ , ^f k _ i iLk 

(2.26)  — - t = -3 — r -f-r  ' 

&V  j -jf  ^V1 


i,k  = l,...,p 


so  that  (2.19)  might  be  written 

(2.27)  ^ 

<bvj<)vK  ^ v j &vk  6vk  6vj 


where  0 « 0 . By  cross-differentiation  as  before  it  is  clear  that 


(2.28)  ^ ^fl  - ^fk  ^fl  * kll 

<)vj  S7  ^ V j l'7k  6-Vj 


for  i j k distinct,  that  is 


(2.29) 


cJv\j  i)1 


,fl  = ^ ^ef  + ^fj  ^ef 
Wj  <5  vk  (Tv^  ^ vj 


Define  a new  differential  operator  on  any  suitably  differentiable 
function 't'  by  means  of 


(2.30)  ^ 


_-V 

’ w 


.... 

1’  ’ n 


h V1 , . . . , ^ -v  n 


I 


! 


Not*  that  although  /'P. 


i i / (^i  - ) 

11>  * ’ ‘ ?1n*l  l,‘",Jn  i 


we  have 


(2.31)  . 


1l’-,-’1n*l 


e^i  i ) 

1*  ’ n ^ 


i i 

n+1 


where  the  customary  partial  derivative  is  indicated  by  a subscript 

without  underlining.  In  particular  ='1*. , /t\.  = 'It.  + 

3 J JK  j K 

and  (2.28)  becomes 

(2 . 321  f1.  - «*P  f1  ♦ fk.  ft  . 

jk  k j j k 


Lemma  C^:  Suppose  'f  satisfies 


(2.33) 


and  let  i^,...,in  be  any  distinct  indices  among  l,...,p,  say 

l,...,n  for  convenience,  n < p.  Then 

n 

(2-3U)  n 


where  1, . . . ,i , . . . ,n  indicates  the  n-1  integers  omitting  i. 
In  particular 

’ k ^ 4 .. 


Proof:  Use  induction  on  n.  By  (2.33)  the  lemma  holds  for  n **  2. 

To  show  that  n may  be  replaced  by  n+1  in  (2.3U)  use  (2.31)  to  obtain 
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1 .... , rv*l 


isf  ) +T  "t 

VT  1. nJ  ^ n*l  'l, n 

—i _ n+1  — 2 <— 


(2.36 


£ • F\,  * T • F1  ♦ 'T  T F1) 

^ v i n*l  i,n*l  n*l  l 7 


n n 

€ 't'iFn.l  * C'+'i.n.l  ?1 


where  F 


1 , . . . ,i , . . . ,n  for  convenience. 


But  the  second  term  on  the  right  is 


U 

<*  t-Ti  fU  *\.i  5*6  f1 


(2.37) 


^*f.  f1  , F1  ♦'Y  ..  f"*1 

<*  l n*l  n+1  1, . . . ,n 


by  (2.33)  and  (2.35),  which  is  just  the  induction  hypothesis  applied 
to  f"41.  Hence  by  (2.3l)  and  the  induction  hypothesis 


15 . . . ,n*l 


= 'Y  (F^  * f^-  F1)  ♦ 'Y  f"41 

% i ' n+1  n*l  ' n+1  1, . . . ,n 


(2.38) 


= i +.  4 a 


+ + £n  X 
1 n*l  1 , . . . ,n 


- £+i  fi  H 


as  asserted 
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Formula  (2.35) 
n directly  in 
If  "+'  satisfies 


may  be  iterated  to  find  an  expression  for 

terms  of  the  f5  if  desired. 

1 

(2.33),  that  is 


(2.39) 


<e+),.  e“  + 

jk 


fj  (e1').."'' 

K 3 


ft  t.1"). 
k 


- 'I' 
e 


then  clearly  0 - e satisfies  (2.27). 
Similarly  (2.3b)  becomes 


(2. bO) 


'V 


,.n 


£ 6 


i=l 


Since  the  right-hand  side  of  (2.bo)  is  symmetric  in  the  indices 
l,...,n  it  is  clear  that  the  computation  of  0.  is  independent 

of  the  order  of  dif fcrentiation.  Thus  any  mixed  derivative  of  0 
involving  no  repeated  differentiation  may  be  uniquely  evaluated  in  terms 
of  the  first  derivatives  of  0. 

THEOREM  C-s  Under  the  hypotheses  of  THEOREM  B,  if  the  integrability 
conditions  (2.23)  are  satisfied,  then  a conservation  law  in  a neighborhood 
of  u - uq  may  be  specified  by  assigning  values  of  0 (=0  ) on  the  p 
characteristic  curves  (2.16)  through  u . In  particular,  by  successively 
determining  conservation  laws  for  which  the  derivatives  of  0 vanish 
along  all  but  the  kLn  characteristic  curve,  k = l,...,p,  it  is  clear 
that  any  system  satisfying  (2.3?)  is  equivalent  to  infinitely  many 
sys terns  of  conservation  laws. 


r 


j £ 

Proof;  Y/rite  0 for  0 (0,...,C),  0 ^ for  0 (0,...,  v , . . . ,0)  , 

0 ^k  for  0 (0, . . . , v, . . . ,v  k, . . . ,0)  etc  . , and  0^  for  > 

>2.  2 k 4v 

0m  for  — j- — , etc.,  and  note  the  following  sequence  of  integral 
fev  & v 

formulae,  which  may  be  verified  simply  by  evaluating  the  integrals: 


(2. Hi)  0^  = 0 ♦ (j  0^  dv^ 


^ k = ( eS  ^ ♦ dK)  - 0 ♦ 


k , 1 


(2.h2)  0^k  = (0  0K)  -*  0 4 j j 0 £k  dv  kdv 

0 0 

(2.U3)  0 = (0kJ  ♦ 0^  ♦ 0 ^k)  - (0^  ♦ 0k  ♦ 0^  ) 


Jvk  J 

04  J J J ^ Jdvkdv^, 


0 0 0 


and  so  forth,  until  0132,...,p  = 0 (V^...,  vP)  at  any  point 


is  piven  bv 


(2.UU) 


0J-5**’*P  = ^(all  0's  with  p-1  entries) 

- (all  0's  with  p-2  entries)  +••* 

, V1  vP 


....  (-1)^0  • J ...  ^ 01; ...f| 


D 1 

dv  ‘ . . .dv 


0 0 


The  data  0^,  i=  l,...,p  is  given,  so  that  ( 2 . Ul)  is  trivial  for 
all  JL  . Everything  on  the  right-hand  side  of  (2.U2)  except  the  integrand 
is  therefore  known.  But  according  to  the  lemma  the  integrand  is  just 


^ 


-38- 


a known  linear  combination  cf  the  first  derivatives  of  0 which  may 

be  integra'  ad  by  parts  to  give  a linear  integral  equation  in  0 for 

which  several  existence  and  uniqueness  proofs  can  be  constructed.  Thus 

0 k 

all  the  0 ‘ can  be  found,  so  that  everything  on  the  right-hand  side 
of  (2.1*3)  except  the  integrand  is  known.  Again  the  lemma  allows  the 
integrand  to  be  uniquely  expressed  in  terms  of  the  first  derivatives 
of  0 which  gives  rise  to  a new  integral  equation,  and  so  forth.  This 
process  successively  finds  the  values  of  0 (V^, . . . ,V^)  on  higher 
and  higher  dimensional  webs  spanned  by  the  characteristics,  finally 
giving  the  values  on  the  p-dimensior.al  web,  that  is,  in  the  entire  space, 
at  least  in  a neighborhood  of  v = vq,  depending  on  the  method  used  to 
solve  the  integral  equations. 
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CHAPTER  III 


. * 


THE  GENERAL  CASE 


§ 1,  Exterior  differential  forms 

In  Chapter  I the  existence  of  conservation  laws  was  entirely 
reduced  to  the  study  of  the  linear  system  (1.23).  This  system  is 
repeated  here  for  convenience  as 

l J \ 0 

(3.1)  J 1 ♦ l2  Uj  02  - 0,  i»l,...,P, 

l 

are  functisns  of  U,  and  the  U ^ are  linear  differential 
first  order  operators,  differentiations  in  the  characteristic  directions 
in  U. 

(3.2) 

* *-  6u 


where 


spanning  the  tangent  space  of  U,  which  is  defined  in  Chapter  I,  ^2. 

The  dual  to  the  tangent  space  may  be  identified  with  the  space  of 
first  erder  differential  forms,  called  Pfaf fian  forms , over  U,  since 
the  former  transforms  ccntravariantly  and  the  latter  covariantly  under- 
a change  of  coordinates,  giving  rise  to  appropriate  bilinear  functionals. 
It  will  be  more  convenient  in  the  present  chapter  to  work  with  the  space 
of  differential  forms,  spanned  by  the  dual  basis 


(3.3) 


u J 


L?l 


dU 


D <j—.  P c" 

where  ($£_)  is  the  inverse  cf  ),  X j 

This  space  is  called  the  co-tangent  space  of  U. 


- £ 


7 

v# 

i 
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An  exterior  differential  form  of  degree  r in  the  variables 
1 N 

z ,...,7  is  any  element  of  the  Grassman  algebra  of  degree  r 

generated  by  the  space  of  Pfaffian  forms  in  the  N variables 
1 N 

z ,...,z  . Thus,  if  the  FTaffian  forms  are  spanned  by  9^,..., 9^, 

an  exterior  differential  form  of  degree  r is  a linear  combination 

with  analytic  coefficients  of  terms  9.  /v  9.  A ...  * 9.  , where 

X1  X2  xr 

i,  are  integers  from  1,...,N,  the  sign  n/^"  denoting 

exterior  product.  The  definition  of  exterior  product  assumes  that 


(3. a) 


0.  A ©• 
*1  2 


9.  A 9. 
12  X1 


= 0, 


These  forms  constitute  a ring  in  which  products  are  formed  by 

exterior  multiplication.  Since  (3.U)  implies  that  the  square  of  any 

element  vanishes,  it  follows  that  this  ring  is  of  dimension 
r N N N 

( _<  ) * ' 0 ) + •••♦(,.)  - 2"  - 1 considered  as  a vector  space 

.1  ci  In 

N 

over  the  analytic,  functions,  thus  of  dimension  2 when  the  unit  element 
is  added,  introducing  the  analytic  functions  as  differential  forms  of 
degree  zero . 

From  ( 3 . U ) one  easily  shows  that  the  forms  G^,...,©^  are 
linearly  independent  if  and  only  if  SiAG2A***A0r  xs  a non~zero 
form. 

The  operation  of  exterior  differentiation,  denoted  by  d,  is 

X1  . xr 

defined  on  a monomial  a dz  a • • • a az 


(3.5) 


X1  1r 
d ( adz  " 


N 

j=l 


^ z j 


']A^ 


p 

! I 


~Ui- 


and  on  arbitrary  exterior  differential  forms  by  means  of 
d(  @ * © ) = d@  + d (©  . Clearly 


i i n 

dd(adz  ^...Adz  r ) = ^ — 


a , k , j , 1 

^ . IT—  dz  Adz  Adz  A- 

j,k<L^Z 


' * A dz 


(3.6) 


= 0 


by  (3.U),  which  gives  Poincare  1 s theorem, 


(3.7)  dd  & = 0 , 

for  an  arbitrary  exterior  differential  form  © . If  & is  an  exterior 
differential  form  of  degree  r and  © is  an  arbitrary  exterior 
differential  form,  then  the  definition  of  exterior  differentiation 
implies 

(3.8)  d((9A©'  ) = d©A®  ' ♦ ( -l)r$Ad  Q ' . 


A differential  form  Cl  which  is  the  exterior  derivative  of  another 
differential  form,  n = d Q , is  called  an  integral,  or  total  differential. 
Poincare's  theorem  has  a converse , namely,  if  d n - 0 in  a neighborhood 
of  a-  given  point,  then  is  an  integral  there,  provided  that  the 
coefficients  of  n are  analytic  in  the  neighborhood. 

The  exterior  derivative  of  any  Pfaffian  form  is  an  exterior 
differential  form  of  degree  two.  Consider  the  forms  uj^,  JL  *■  l,...,p 
given  in  (3.3)  for  example.  We  have 


dw 


i 


U J 


(3.9) 


where  we  rai^ht  as  well  take  c^g-  ♦ c^p  = 0 due  to  (3.U). 
If  f is  an  arbitrary  analytic  function  of  u\...,u^  then 


(3.10)  df 


if  , <r  y^if  -tP,  x 
- du  = (a  Jydu 


V f UJ 


so  that 


0 = ddf  = dCU^  f)AuJ°"+  Urf  du or 

- (U,f)/AA  Ux  f di^r 


(3.11) 


" \ {u  f>  (Ug-  f)  - IV  (Up  Ur  f d^ 

= i [UpU^lf^^  Urf  d^1 

l„t  ? r 1 ~c  P <r 

^•itfptrU-  f UJ  alO  _ | CpV  ur  fuJ> 


•p  ( ^p<r  - c p o-  ) U^  f uJ  */4  u.> 


P .,>0' 


by  the  definitions  (1.29)  and  (3.9)  of  the  structure  functions 

yt  - 

flp°"  and  the  functions  CpV-  respectively.  3y  definition  opa-  = c 
for  P = 0‘  ; and  since  (3.11)  implies  that  the  coefficient  of  ’ ^ ^ caj tT” 
must  vanish  for  O’  7^  p it  follows  that  the  Cp<r  of  (3.9)  are  merely 
the  structure  functions. 


(3.12) 


- icXT 


,r 

pv  , 


f.CT  ,T  = l,...,p. 


I 


§ 2 . Frobenius 1 theorem 

Tn  Chapter  II,  Jl,  it  was  noted  that  the  complete  system  (2.2), 


repeated  here  as 


(3.13) 


- 0, 


= 1 , . . . , l 


for  convenience , where  W n * ^ f jj  — — is  a linear  homogeneous 

^ £>wi 

1=1 

first  order  differential  operator  in  a neighborhood  of  w^  in  a 
manifold  W of  dimension  N,  possesses  N-r  functionally  independent 
solutions.  It  was  further  noted  that  there  is  a change  of  variables, 
w = w(v) , such  that  the  given  system  becomes  Jacobi  complete, 


(3.lU) 


6$  _ n 


H. c j 


The  solutions  (£>  of  this  system  are  arbitrary  functions  of  v ^ , 

^ = r+1 , . „ . ,N.  Thus  one  can  map  a neighborhood  of  itQ,  given  by  the 
coordinates  vq  = (v^,...,v^),  onto  a neighborhood  of  some  point 

m = (mrt^...m^)  in  an  (N-r) -dimensional  manifold  M by  means  of 
o o’  ? o 


(3.15) 


J - V i 


r+1, . . . ,N. 


M is  the  integral  manifold  determined  by  1,3.13)  in  a neighborhood 


of  w 


The  fact  that  the  vectors 


r+1, . . . ,N  span  the  tangent 


space  of  M in  a neighborhood  of  is  of  no  interest  since  these 

vectors  depend  on  an  a priori  knowledge  of  M.  However,  the  differentials 
dv^  = dm1^  / = r+l,...,N  span  the  corresponding  co-tangent  space  and 


-kh~ 

can  be  invariantly  described  as  isomorphic  bo  the  annihilator  of  the 
vectors  , -I  - l,...,r,  given  in  (3.13)  and  (3.lU),  which  determine 

subspaces  of  the  tangent  spaces  in  a neighborhood  of  w . Furthermore, 

since  these  differentials  are  independent  their  exterior  product  is 

non-zero,  and  so  one  could  equally  well  describe  M as  the  integral 

r*l  N / 

manifold  through  w for  which  dv  A ...Adv  f 0. 

It  is  desirable  to  find  a criterion  for  the  completeness  of  (3.13) 

in  terms  of  the  annihilator  A of  Suppose  (3.13)  is 

complete,  and  let  ©,, , « • • ,0^  be  any  basis  of  A.  Then  each  0 £ is 

r+1  N r*l  N 

some  linear  combination  of  dv  , ...,dv  , where  v , ...,v  are  any 
independent  solutions  of  (3.13),  with  coefficients  which  are  not  a 
priori  known.  Since  the  square  of  any  Pfaffian  form  vanishes,  only  one 
term  appears  in  the  product  -O-  = Q ^ A • • • aSN’  name-*-y> 

(3.16)  n-  ©r+1  A ...  a0n  - A dvr+1  A ...  A dvN. 

Thus  A is  completely  described  by  a non-zero  form  i"!  for  some  unknown 

factor  /\  , A t 0”  Differentiating  the  right-hand  side  o.f  (3.16) 

1 r+1  N 

according  to  (3.6),  all  cerms  except  dA^dv  A . ..^dv  vanish 
by  Poincare's  theorem.  Hence 

(3.1?)  d_fL  = ^ a -I~l 

so  that  d _Tl  is  in  the  ideal  generated  by  J\  over  the  ring  of  differential 
forms,  which  is  also  expressed  by  saying  that  the  differential  form  -T\- 
is  closed.  Conversely,  if  -f).  is  closed,  Frobenius 1 theorem  asserts  that 
there  is  an  integral  manifold  through  for  which  -Tl / 0,  that  is, 

satisfying  (3.13). 


t 
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As an  application  we  return  to  the  problem  of  Chapter  IT,  J1  to 

determine  when  multiplication  by  an  integrating  factor  and  a change  of 

_ into  the 
iu^ 


variables  u ■»  u(v)  will  bring  the  operators  U«  5^, 


form  A ^ s 7^1  . If 


(3.18) 


7 


then  the  dual  basis  is  of  the  form 


(3.19)  ou^  2.  du^  =A^dv^,  /*»  l,...,p 

0 JL 

Thus  for  each  -t.  the  differential -H.  = wJ  must  be  closed.  That  is, 
the  exterior  fonn 

(3.20)  duo1  - - | 

p J. 

must  lie  in  the  ideal  generated  by  u>  , so  that  must  vanish  except 
when  o~  ^ A or  . This  condition,  found  previously  ir  Chapter  II, 

is  both  necessary  and  sufficient  by  Frobenius ' theorem. 


3.  Exterior  differential  systems 

Instead  of  prescribing  only  'a  form  -TL  which  does  not  vanish  on  some 
unknown  integral  manifold  one  might  also  seek  integral  manifolds  on  which) 
in  addition,  a prescribed  system  of  exterior  differential  forms  does  vanish. 
Problems  of  this  type  arise  in  a natural  way  from  systems  of  linear 
homogeneous  first  order  partial  differential  equations.  For  example,  it 
is  of  interest  to  consider  the  system  of  Chapter  I,  £5,  concerning  the 
case  p = 2, 


■21)  ^ i o 4 $2  ri  = 1)2 


/t.  | 2 

where  T . are  known  functions  of  v such  that  „ _ 

J 1 }•  2 r 2 


si  sj 

^22^  °* 


For  simplicity  we  introduce  the  unsymmetry  of  Chapter  II  and  take  this 
in  the  form 


(3.22)  ^ =0,  i-1,2. 


<)v^  iv^ 


(3.23)  d02  - 0-,  dv1  4 02  dv2 


then  (3.22)  implies 


(3.2U)  d0d  - - 5 10^dv1  - j202dv2. 


Suppose 


(3.25) 


d0x  *=  01]L  dv1  ♦ 012  dv2 


d02  = 021  dv1  4 022  dv2  , 


then  exterior  differentiation  of  (3.23)  gives 


0 = (011dv1  ♦ 012dv2)_A  dv1  4 (021dv1  4 022dv2)  A dv2 


(3.26) 


012dv2A  dv1  4 021dv1Adv2 


22  'a 


m (021  ~ ^12}  dylA  dy2 


L 
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so  that 


(3.27)  021  - l\2 ° 

u * h.  20  - ^ 

This  tos  to  be  expected  since  it  merely  states  that  Xv2  Av1 

Similarly  exterior  differentiation  of  (3.2U)  gives 

(3.28)  5 j02  * 3 2021  - I X * 5V12 


where 


.31  i^-dv1  ♦ 3V,  i - 1^2 


Hence,  since  by  assumption 


I2  i 31, 


(3.29)  021  * 0i2  * AX  + A 02 


1 J2 


2 51 


where  /J-  - x2  ^1  and  A y 


Now  substitute 


these  values  into  (3.25)  and  take  exterior  derivatives  to  find 


(3.30) 


0 w d0-,  Adyl  4 P^Adv2 


0 « |°2dv2Adv1  ♦ d022- Adv 


where 


(3.31) 


f 1 = yW-^u.  4 CAi  ♦aV  )01  + Xl  + A2  ,02 


p2  - /022  + (/U2  4/AV  )01  * (A  4 aV  )02‘ 


Finally  (3.23)  and  (3.30)  may  be  taken  together  as  a differential 
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system  in  the  four  unknowns  0^,  02,  $22’  and  two  indePenden't 

variables  vd  and  v2 , 

/ d01  - 0ndv1  - (/u.'L01  + A202  )dy2  ' 0 
d02  - (y/.^  ♦ ^$2  ^dyl  ” ^22dy2  ■ 0 

(3.32)  J 

d0i;L  A dv1  ♦ |°  1dv1A  dv2  = 0 
\ P ^ dv2A  dvd  ♦ d022 a dv2  “ 0 . 

The  relations  d(0^dvd  ♦ 02dv2)  = 0 and  d ( J 'L01dvd  ♦ J 202dv2)  * 0 

imply  that  the  two  Pfaffian  forms  in  parentheses  are  integrals,  by  the 

converse  of  Poincare's  theorem.  Since  these  two  equations  can  be  obtained 

without  differentiation  from  (3.32),  we  have  omitted  (J.23)  and  (3.2U) 

from  the  collection  (3.32),  In  fact  this  system  already  describes  the 

1 2 

original  unknowns  0 and  0"  up  to  additive  constants  of  no  interest 
in  conservation  laws. 

Now  note  that  further  exterior  differentiation  of  (3.32)  yields 

only  equations  which  are  already  in  the  ideal  it  generates  over  the 

1 2 

ring  of  exterior  forms  in  the  space  (v  ,v  *0qj02,^ll,^22^ ’ so  that 

the  system  (3.32)  is  closed.  If  an  integral  manifold  can  be  found 

1 2 

which  satisfies  (3.32)  and  on  which  in  addition  dv  ^ dv  0 0,  then 

the  system  is  said  to  be  in  involution  with  respect  to  the  variables 

12'  1 
v and  v . Clearly  if  (3.3 2)  is  in  involution  with  respect  to  v 

O 

and  v then  there  exist  solutions  of  (3.21).  The  number  of  such 


* **  f V*  rr  's~-^ 


i 
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solutions  as  well  as  a criterion  for  a closed  system  to  be  in  involution 
is  given  in  §h. 

An  important  detail  in  the  preceding  derivation  should  be  stressed. 
The  original  differential  system,  comprising  (3.23)  and  (3.2h),  was 
not  closed,  but  it  might  have  been  closed  merely  by  adding  the  exterior 
derivatives  of  these  equations  to  it,  giving  the  system 


(3.33) 


/~  a02  - 0,dv1  - 0odv2«  0 

/ 1 2 

d0^  * !"''0.]dv'*'  * } 20gdvc  * 0 

) 

\ 1 2 

d0^  A dv  ♦ d02  ^ dv  **  0 

\ j 1¥1  a dvl  ♦ I 2d02  a dv2  ♦ (02  } 2 - 0152  )dv1  A dv2  - 0, 


in  which  the  first  two  equations  might  be  omitted  as  before.  As  it  turns 
out,  (3.33)  is  not  in  involution,  so  that  it  has  been  necessary  to  prolong 
it  by  adding  0.,  and  02o  to  it  by  means  of  (3.25).  Then  the  new 
system  had  to  be  closed  by  the  addition  of  the  exterior  derivatives  of 
(3.25) J luckily  the  resulting  system  is  in  involution  as  we  shall  see  later 
on.  It  is  an  open  question  whether  an  arbitrary  system  can  be  prolonged 
in  this  fashion  into  a system  in  involution  with  respect  to  a given  set  of 


variables . 

The  remainder  of  this  section  is  devoted  to  replacing  the  general 
equations  (3.1)  governing  the  existence  of  conservation  laws  by  a- 
closed  exterior  differential  system.  Since  one  may  solve  (3.1)  for 
01  in  terms  of  some  new  parameters,  say, 


r 


(3.3U) 


"f  ^ ■ XIJ2 


U if  " , 


the  initial  differential  system  in  the  variables  (0^,0^,X-^,...,X  ,u\...,uP) 


(3.35) 


d01  - ^ X p ^ “ 0 


,2  < i 

d0  ♦ f X,31  ^ - 0 


. P j 

where  = . du“.  The  system  (3.35)  can  also  be  obtained  directly 

> J 
j=l 

from  (1.9)  without  going  to  the  trouble  of  deriving  the  corresponding 
linear  equations  (3.1).  In  fact,  write  (1.9)  in  the  normal  form 

*£  i ,i  i n 

(1.18)  j 2 °”i  “ -3  1 °~2  = °»  X«l,...,p 


l ^ l hx? 

where  0~  . = S \ • - — - , and  suppose  there  exists  a linear 

1 r;  J ^xx 

,i*i 


combination 


y H 0 £ 

^X.  ^ (J  3 £T-^  - I i°^2  ^ -*-e^~hanc*  members  of 


(1.18)  which  is  of  the  form  £ *K-  Then 

dx  &x 


dx^  l 


°~1 


<>  X*  * 

2 - ^X/Sl°V 


since  each  of  these  identities  involves  differentiation  with  respect  to 


only  one  of  the  independent  variables,  we  may  replace  the  derivatives 

and  <T  f by  the  corresponding  differentials  d0^  and  w ^ to  obtain 
6 x1  1 

(3.35) . 


We  attempt  to  prolong  (3.35)  by  adding  new  variables  to  it, 

where 

p 

(3.36)  dX.^  =*-j2pUJ  > i=  1,...,P, 

Denote  the  Pfaffian  form  on  the  right-hand  side  of  (3.36)  by  Oj^  . 

The  introduction  of  0 ^ is  convenient  for  deriving  results  which  are 

independent  of  the  fact  that  it  is  an  integral;  it  should  be  noted  in 

/ 

particular  that  Poincare's  theorem,  dd  Qj  = 0,  is  always  valid,  even 
if  nothing  is  known  about  QI- 

Differentiate  (3.35)  to  obtain 


■ 3.37)  ^ { QJlAf  d ($*  w*)j  - 0,  i - 1,2. 


The  coefficient  of  in  (3.37)  is  given  in  terns  of  the  X» 

H x % 

and  the  known  functions  in  the  second  term,  and  by 

X d J k - X„,  i . in  the  first  term,  i *>  1,2.  Hence,  since  the 
k * i i 

characteristic  directions  in  X are  distinct,  by  setting  the  coefficients 

0 If  A 

of  u>  UJ  equal  to  zero,  k,  X-  = l,...,p,  one  can  solve  for  all  those 
X j|  for  which  k / Jl  in  terms  of  the  %£  . Taking  these  values  for 
y.  k ^ , or,  preferably,  merely  considering  these  values  as  definitions  , 
(3.37)  is  identically  satisfied  and  hence  the  exterior  derivative 

( 3 - 38)  jdO^l^-  9iAd(^uJ£)  ♦ dXjA  £0 

i-  i 1 


! 
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of  the  left-hand  side  is  identically  zero.  (Note  that  a minus  sign 
appears  in  (3.38)  due  to  (3.8).)  In  other  words,  since  the  la3t 
two  terms  of  (3.38)  cancel,  using  (3.36)  and  the  definition  of 
Q , the  expressions 

p 

(3.39)  S dO„A}*u^ 

A ^ 

& =1 

are  in  the  ideal  generated  by 


(3-Uo)  dXn- 


= p. 


For  the  converse  of  Pbincare's  theorem  turns  (3.35)  into  (3.37), 
and  (3.37)  vanishes  identically,  merely  being  used  to  define  various 
coefficients  in  the  Qj's.  Note  that  the  system  (3.U0)  is  not  yet 
closed,  although  9^  represents  an  integral,  since  © ^ involves  new 
terms  which  prolong  the  original  system.  In  the  prolonged  system  0 j 
is  not  an  integral. 

Now  let 

( 3.Ul)  dQ  ^ - Gj^auA  | 

for  some  , TT^V^  depending  on)u|  ^jj  and  their  first  differentials, 

/,  j,k  = l,...,p,  where  TT^  ♦ = 7et  closing  (3.U0)  by 

Q 

asserting  t’nat  dO  ^ - 0,  multiply  (3.Ul)  by  to  find  that 


(3.U2)  J £ j • 

4j>k=i 

lies  in  the  ideal  generated  by  (3.U0)-,  since  (3.39)  does; 


i 

f 

i 
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hence 


(3.U3) 


•*i  ^.Ik  ^jkl  ^i^kij  " -°»  1 = 1>?' 


But  (3.^3)  gives  the  identity 


(3-1*1*)  TT; 


si  Jj 

5 » 55 


JLisZ 

Tk  ri 

3 1 5 1 

Tk  ^ 

i 2 3 2 


TT(ijk) 


J 1 ; 1 


I r-^  rJ 
3 2 32 


where  IT  ( £ jk)  is  independent  of  the  order  of  ^ jk.  In  the  special  case 
of  Criapter  II  this  simple  result  is  exactly  the  assertion  of  Chaptei  II, 

^ 2,  that  only  one  first  order  relation  arises  for  each  distinct  triple 
of  indices.  In  fact,  TT(i jk)  is  just  the  left-hand  side  of  (2.37).  It 
should  be  noted,  however,  that  in  general  IT  is  a linear  combination  of  the 
'X-Jl  as  well  as  X|  , Jl-  l,...,p,  the  special  case  being  the  only 
exception  to  this  rule. 

It  is  undoubtedly  true  that  for  p > 3 one  can  find  still  further 
identities  among  the  If  ( ^jk).  These  identities,  if  any,  will  not  be 
investigated  here. 


The  fact  that  dd Qj  belongs  to  the  ideal  generated  by  (3.1*0), 


that  is 


(3.1*5) 


d ? 


I A ir»,ic 


j,k  t* 


2 <1  fjk  d ^ ~ °> 

J jk,  // 


will  be  recorded  here  i'or  later  use* 


To  close  (3.U0)  one  must  first  add  to  it  the  requirement 
d9^  * 0,  which  breaks  into  several  parts  in  the  notation  of  (3.Ui) 
The  result  is  that 


(3.U6) 


TT  ( £ , j , fc)  ' 0 

~ ®i  " 0 

^ 0 

OU,)  4 A 0»J  * =0 


/, j,k  - 


is  a closed  system  provided  that  d(uJ^  ) and  dT"  ( H jk)  belong 
to  the  ideal  generated  by  it.  However,  it,  is  not  necessary  to  look  at 
d(uh^AU>^)  since 


d(u3^Auu^)  ♦ \ d1T  /jk  u'^Aa,k> 

JS-M 


belongs  to  the  ideal  generated  by  the  left-hand  members  of  (3.^6), 
thanks  to  (3. 45).  Hence  d(u)^AU^)  vanishes  if  dTnijk)  does 
for  all  j , k — l,»»»,p. 

Finally,  if  dTT(^jk)  vanishes  then  so  does  its  derivative,  by 
Fhincare's  theorem,  since  no  prolongation  of  the  system  is  involved  at 
this  point.  Thus 


(3.U?) 


• 7T  ( i , j,k) 

= 0 

dTT{J  , j,k) 

“ 0 

i.  * 1,...,P 

dlLl  ~ QJt 

= 0 

f 

= 0 

L__ 
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is  a closed  system  in  X^  ,X^,  and  u^  , .i  * l,...,p,  all  of  whose 
solutions  are  also  solutions  of  (3.1).  The  main  result  of  this 
chapter  is  that  (3.U,7)  is  also  always  in  involution  with  respect  to 
the  u ^ . 


L 


§U.  The  theorem  of  Cartan  and  Kahler 

The  following  existence  theorem,  due  to  Cartan  [ 1 ] and  extended 
by  Kahler  [2  ] t.o  arbitrary  exterior  differential  systems,  not  only 
gives  a criterion  for  a closed  system  to  be  in  involution  at  a given 
point,  but  also  indicates  how  much  Cauchy  data  (in  the  form  of 
coefficients  of  convergent  power  series  in  one  or  more  variables  about 
the  point)  is  necessary  to  determine  a unique  solution  of  the  system. 
The  presentation  given  here  for  the  case  that  the  system  contains  no 
forms  of  degree  higher  than  two  is  essentially  that  given  in  Kahler, 
pp.  5^-55,  with  only  slight  changes  of  notation.  Frobenius'  theorem 
and  the  converse  of  Poincare's  theorem  can  be  obtained  as  special  cases 
of  this  theorem. 

Given  p independent  Pfaffian  forms  uA,...,luP  in  a certain 
neighborhood  in  some  space  of  dimension  n * r,  we  consider  a system 


TT.  - 0, 


(3.U8)  / 9.  - °, 


l Xk-0, 


i s 1 j • • • 


j » l,...,h 


k 1, . . . ,m 


of  forms  TT.,0  X,  of  degree  zero,  one,  and  two  respectively,  where 
l j k 
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the  Jacobian  matrix  formed  by  the  gradients  of  the  forms  of  degree  zero 
is  of  rank  r,  the  p * h Pfaffian  forms  ©^,...,©^  are 

linearly  independent,  the  quadratic  forms  are  of  the  type 
P . 

x*~l*^*  for  some  new  Pfaffian  forms  uj  ^ , and  the  coefficients 

i«i 

of  all  of  the  forms  or“  analytic.  *Ve  seek  a solution  of  (3.U8)  on 
which  K / 0. 

Let  q = n - p - h and  suppose  that  new  independent.  Pfaffian 

1 P 

forms  are  chosen  in  such  a way  that  wJ  0 , 6^,...,©^, 

span  the  co-tangent  space  of  the  n -dimensional  manifold 
determined  by  the  zero-dogres  forms 


(3.U9) 


“HV  ”0,  i = 1,... ,r. 


Then  the  Pfaffian  forms  are  linear  combinations  of 

1x1  = (tu  , . . . ,u>  **)  , © = (Q^ , . . . ,9^)  and  w **  (Co^, . . . , Cu^);  s ince 
only  the  coefficients  of  uj  are  of  interest,  and  only  the  manifold 
determined  by  setting  TT  = (TT  ^....IT  ) = 0,  we  may  write 


(3.50) 


* 2 °iis3 


. (mod  u;,©,  TT  ) 


k = 1 , . . . ,m;  ^ ” l....,p. 


Let  represent  the  number  of  linearly  independent  forms  among 


Co  ^ u^,  mod  uJ  ,0,  TT  , k = 1, . . . ,m,  where  u^  * (u^, . . . ,u^) 

4-i 
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a set  of  constants  chosen  in  such  a way  that  0"%  is  a maximum , or, 

preferably,  u^  is  a fixed  set  of  indeterminates . Similarly  let 

crn  ♦ CT„  represent  the  num.ber  of  linearly  independent  forms  among 
Xc 


Hju!  and  i 

i -i  i “i 


0 

^ u2  mod  UJ  ,©,  Tf  , 


where  u^  is  chosen  in  such 


a way  that  Is  a maximum,..., 


let  rr  +...♦  <r 

1 P-1 


represent  the 


number  of  linearly  independent  forms  among 


l 

i “i 


mod  ^ ,6.  TT  , 


where 


u 


P-1 


is  chosen  in  such  a way  that 


0"  , is  a maximum. 

P-1 


Finally  let  CT 

P 


q 


<rp-l* 


One  further  definition  is  necessary.  An  integral  of  (3.b8)  on 
which  ...Au>  f 0 is  completely  specified  once  w is  written  in 

terms  of  uJ  , say, 


(3.5D 


essentially  prolonging  (3.^8)  by  adding  the  t^ 
t^  are  restricted  by  the  quadratic  relations  in 
in  general  span  a space  of  pq  dimensions.  Let  M 
the  space  spanned  by  the  t ^ . 


to  it.  However,  the 
(3.U8)  and  so  do  not 
be  the  dimension  of 
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Ths  Cartan-Kahler  theorem  asserts  that 


(3.52)  M<  pq  - (p-l)«r1  - (p-2)<r2-...-<Tp_1 


with  equality  holding  if  and  only  if  (3.^8)  is  in  involution  with 
respect  to  p variables  for  which  u>^ A ...  f 0.  In  the  latter 
the  general  solution  depends  on  cr  arbitrary  parameters,  where 


C~  = n - p - q,  on  0“  arbitrary  analytic  functions  of  a single  variable, 

O 1 

on  (T~2  arbitrary  analytic  functions  of  two  variables,...,  and  on  CT 
arbitrary  analytic  functions  of  p variables. 

This  theorem  is  proved  by  means  of  a sequence  of  Cauchy- Kowalewski 
constructions.  Setting 


uO- 


~ u^u;*  one  successively  constructs  "integral  elements* 

k=l 


through  the  spaces  determined  by  1 =...*  = 0 for  k = l,...,p. 

It  should  be  remarked  that  even  though  a closed  system  is  in 
involution,  further  prolongation  might  still  lead  to  a new  system  in 
involution.  In  particular,  the  singular  solutions  of  a system  in 
involution,  which  satisfy  additional  restrictions  under  which  the  given 
system  is  no  longer  in  involution,  can  sometimes  be  obtained  as  the 
general  s olutions  of  a new  prolonged  system  in  involution. 


ine  existence  of  conservation  laws 
As  a simple  illustration  the  Cartan-Kahler  theorem  can  now  easily  be  applied 
to  the  first  example  of  £ 3,  given  by 
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(3.32) 


'd01  - 011dv1  - (/tV-L  ♦ /M-202)dv2 
d02  - {jj.  101  /u202)dv1  - 022dv2 

(d0H  - f xdv2)A  dv1  = 0 


0 

0 


(d0 


22 


- p2dv1)Adv" 


0. 


Since  there  are  no  forms  of  degree  zero,  r *=  0.  The  variables  occuri-ing 
1 2 

in  (3.32)  are  v , v , 0^,  0?,  0. , , 022  so  that  n = 6.  A solution  is 
1 2 / 

sought  on  which  dv  ^dv  f 0,  so  p = 2.  The  system  (3.32)  contains 
two  independent  Pfaffian  forms,  giving  h = 2 and  q=*n-p-h*2. 

Clearly  we  may  set  = d0..j,  with  = UJg  (mod  uJ,Q),  X = 1,2; 

thus  0“^  * 2 and  0"2  ■=  q - * 0.  Finally  d0^  and  d022  are 

12  1 2 
completely  determined  by  d0-^  * t-^dv  * /°]_dv  and  d^22  * (°2dv  * 't22<iv  5 

where  t-^  and  t22  are  unrestricted  by  the  quadratic  members  of  (3.32), 

giving  M = 2.  Thus  equality  holds  in  (3.52)  and  the  general s olution 

of  (3.32)  depends  on  two  parameters,  and  two  functions  of  a single  variable, 

since  n - p - q * 2 and  <T“  •>  2. 

More  generally,  the  special  case  of  Chapter  II  follows  exactly  the 

same  pattern.  Again  there  are  no  forms  of  degree  zero,  - 0,  provided 

the  integrability  conditions  of  THEOREM  B are  satisfied.  Referring  to 

(3.U7),  it  is  clear  that  each  contains  only  one  term  of  interest, 

namely  d>.»^  , which  was  defined  in  (3.36),  so  tfiat  we  may  take 

with  uj •jq  = uJ£  (mod  u)  ,9) , X = l,...,p.  Then  exactly  as 

before,  n = 3p,  h « p,  q * n - p - h = p,  (T^  = P,  <T2  ». . .*  0^  - 0, 

since  CT  ^ <J~  = q,  and  M = p.  The  equality  holds  in  (3.52) 


H 


•T* 


J * 
i 

i 

I 

f t; 

» * 
i * 
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and  so  the  general  solution  depends  on  p parameters  and  p functions 
of  a single  variable.  This  result  might  be  slightly  surprising  since 
the  existence  theorem  given  in  Chapter  II  indicates  that  the  values  of 
the  first  derivatives  of  X along  the  p characteristic  curves  determine 
a solution  uniquely  up  to  two  additive  constants  which  do  not  occur  in 
the  present  discussion.  However,  if  one  assigns  not  the  first  but  the 
second  derivatives  along  the  characteristic  curves,  the  p values 

of  the  first  derivatives  at  the  origin  may  be  taken  as  the  additional 
parameters . 

As  a slightly  more  interesting  example  consider  the  special  case 
again  but  suppose  that  none  of  the  integrability  conditions  is  satisfied. 
Then,  since  the  coefficient  of  occurring  in  IT  ( ^jk)  is  non-zero, 
the  equations  dTT  jk)  =0  in  ( 3 - U? ) provide  a determination  of  all 
the  in  terms  of  the  , the  determination  being  unique  since  (3.U7) 

is  closed.  Because  the  system  already  contains  expressions  for  d/*j  , 
namely  d'^  = , the  quadratic  members  of  (3.U7)  are  redundant  and 

may  be  dropped  out.  When  no  quadratic  terms  are  present  both  sides-  of 
(3.52)  vanish  so  that  the  system  is  automatically  in  involution.  No 
arbitrary  functions  may  be  assigned,  since  q * 1,  and  the  number  of 
arbitrary  parameters  depends  on  the  equations  TT  ( jk)  = 0.  The  number 
of  these  equations  is  ( ? ),  but  there  might  be  considerable  linear 
dependence  among  them,  as  suggested  on  p.  53.  However,  for  p = 3 
there  is  exactly  one  relation  among  the  three  unknowns  so  that  a 
two-parameter  family  of  conservation  laws  exists. 

To  illustrate  the  preceding  result  consider  the  system 
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which  satisfies  none  of  the  conditions  of  THEOREM  B.  In  fact  it  leads 
to  the  first  order  equation  TT  ■-  o w here 

crM  rr  of  2 3^02  * / 1 2 A 3x  . , 1 2 3,  & 02 

(3,5a)  II  = 2 ( v - v )^T  ♦ ( - v - v + v ) + (v  - v + v ) 

V V ^ V 


Nevertheless  it  has  two  solutions, 

(3.55) 


*1  1,3  2, 

0 = V (v  - V ) 


02  - v1  + v2  + v3' 


and 


(3.56) 


/ wl  1,3  2x  , 1 . 2 A 3v 

( 0 = V (v  - v)(v  ♦ V ♦ V ) 


-,2  , 1.2  . , 2.2  A , 3*2  A 1 - 2 A 3, 

0 = (v  ) + (v  ) ♦ (v  ) ♦ v (v  ♦ v ) , 


the  general  solution  being  a linear  combination  of  these  with  constant 
coefficients . 

Intermediate  cases,  in  which  some  but  not  all  of  the  integrability 
conditions  of  THEOREM  B are  satisfied,  can  be  investigated  in  exactly  the 
same  fashion.  In  every  case,  of  course,  the  conservation  laws  fall  3hort 
of  equivalence  to  the  original  partial  differential  system  unless  all  of 
the  integrability  conditions  are  satisfied. 
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Some  examples  not  contained  in  the  special  case  Tri.ll  be  given  in 

order  to  motivate  the  theorem  to  follow.  In  giving  these  examples  it 

will  be  convenient  to  specify  the  differential  forms  . . . , or 

equivalently,  the  operators  , merely  by  giving  their  structure 

functions,  wh^ch  will  be  taken  as  constants  satisfying  the  quadratic 

relations  arising  from  Jacobi's  identity,  (1.32).  The  fundamental 

1 P 

theorem  of  Lie  groups  insures  the  existence  of  w , . . . , uu  in  an 
appropriate  neighborhood.  Y/e  take  p **  3 and  suppose  that  the  ratios 
^ \ 2 are  distinct  constants.  Further  we  suppose  for  simplicity 

that  the  structure  constants  satisfy  a condition  which  is  just  the  opposite 
of  that  assumed  in  Chapter  II,  namely  that  J is  non-zero  only  when 

JK 

i j k are  all  distinct.  For  convenience  let  (ijk)  be  any  even  permutation 
of  (123)  and  write 


(3.57) 


»]k 


and 

(3.58) 


— 


X . 

i 


U tk 
J 1 >1 


■>2  ->  2 


Then  differentiating  (3.35)  gives 
(3.59)  dXi-01-X  ♦ 


X . *k 


and  further  differentiation  gives 


L 
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(3.60)  Tf*  = iXXi±*  X jj  ♦ ?kXkk 
and 

(3.61)  Au;i  - o,  i - 1,2,3 

where 

(3.62)  to. . = d X . , - 2 tfk  ( y.  ujj  ♦ * uJk). 

' li  ii  j k 


Now  we  consider  various  examples  obtained  by  substituting  these  values 
of  ©.  , UO  . . , and  TT  into  the  closed  system 

i*  ii  ’ 


(3.U7) 


(TT  * 0 
dTT  = o 
d Z±  - = 0 

= 0, 


for  several  choices  of  6 \ tf^,  *1  The  special  case  tf ' •»  if  ^ 


is  omitted  here. 

1 2 

Example  I:  Suppose  tf  =1,  tf 


0.  Then  tt>  is  known  in 


12  3 

terms  of  the  eight  variables  u ,u  ,u  , Xg,  X^,  s^nce 

(3.60)  implies  * 0.  Hence  there  are  only  two  quadratic  relations 

in  (3.6l).  Furthermore  q=n-p-he8-3-3’2  so  that  since 


is  clearly  2,  0”^  = = 0.  Since  the  components  of  ^*^22 

2 3 

and  d"t^  in  the  directions  uj  and  tu  respectively'  may  be  arbitrarily 
prescribed,  M = 2.  Hence  M = pq  - (p-l)0"^  ”2  so  that  the  system  is 
in  involution.  Furthermore,  three  parameters  and  two  analytic  functions 


) 


of  a single  variable  determine  the  corresponding  conservation  laws 
since  n - p - q = 3 and  a~  = 2.  It  is  interesting  to  note  that 
d X ^ * 0,  so  that  X ^ is  one  of  the  parameters.  A more  careful 
examination  of  the  Cartan-Kahler  theorem  would  probably  show  that  the 
values  of  X ~ and  X ^ at  a fixed  point  are  the  remaining  parameters, 
giving  rise  to  an  equivalent  system  of  conservation  laws. 


12  3 

Example  II:  Suppose  If  = & =1,  If  J = 0.  Here  we  note  that 


u-’  = d X ,,  and  = d 'X. 

IX  XX  cd 


22  so  that  dfT  *»  0 implies  ^2  = 


Since  the  quadratic  relations  show  tbat  and  lie  in  different 

X X • dd 


directions  in  the  co-t.angent  space  this  means  oj^  « *22  * ^-eav^-nS 

only  one  quadratic  relation,  uj  ^uP  ” 0»  As  before  there  are  eight 
variables,  which  may  be  taken  as  u ,u  ,u  , X^,  ^22’  ^33 


to  the  relation  TT  =•  0.  In  a fashion  similar  to  the  first  example  one 
computes  M = q - = 1,  where  it  should  be  noted  that  the  additional 

Pfaffian  relation  dTT  =0  gives  h = U rather  than  h = 3.  Thus 
M = pq  - (p-ljff^  = 1,  so  that  the  system  is  in  involution.  One  may 
prescribe  three  arbitrary  parameters  and  one  analytic  functions  of  a 
single  variable  in  order  to  determine  the  solution,  since  n - p - q = 3 
and  <T  = 1. 

Example  III:  Suppose  If  ± - ^ «*  o ■ 1.  In  this  case 

WJ  = d Z ..  - 2(  ♦ X.  W k)  according  to  (3.62) , for 

n n j k 5 b ’ 

i =.1,2,3,  so  that  the  Pfaffian  equation  dTT  = 0 becomes 

(win  ♦ £>22  * ^33  ) * **(  ,Z1wjl  * *3^  ) * o.  Since  LU  ±± 

is  in  the  direction  of  wi1  this  breaks  into  three  Pfaffian  relations. 


Vi# 
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name  ly  u5..  ♦U.w1*0,  i = 1,2,3.  Thus  all  the  UJ . . are  known, 
and  the  quadratic  relations  are  merely  redundant.  Note  that  although 
six  Pfaffian  forms  are  involved  in  the  system  itself  only  five  of  them 
are  independent  since  the  relation  d TT  =0  is  merely  the  derivative 
of  the  definition  of  X--q,  say,  which  is  not  a prolongation.  Thus 
clearly  q = n-  p-  h = 8-  3-  5 = 0 and  M * 0~  ^ = 0 so  that  the 
system  is  in  involution  and  is  completely  determined  by  five  parameters 
since  n - p - q = 5- 


In  every  example  so  far  the  relations  m = M = = q and 

C3  (j—  =0  have  been  satisfied,  where  m is  the  number  of 

independent  quadratic  relations,  so  that  the  condition 

M = pq  ( p>— 1)  CT <y  , has  always  been  satisfied  by  virtue 

1 P“1 

of  m = pm  - (p-l)m.  If  is  clear  that  = m always  holds  since 

each  independent  contains  exactly  one  term  independent  of 

the  Pfaffian  forms  . , uU and  the  other  For  the 

same  reason  C~^  “ ...  = Furthermore  M = m always  holds 

o 

since  the  component  of  dX^  in  every  direction  except  is  determined 

for  each  term  occurring  in  the  m independent  quadratic  relations;  that 


is,  in  the  language  of  ^U,  there  are  n independent  parameters  . 

Finally,  to  show  that  q = m in  every  case,  suppose  the  opposite,  that 

•I 

is,  q > m since  q < m is  impossible  by  the  Cartan-Ka'nler  inequality 


(3.52) . In  this  case,  since  q - 0“^  ♦ ...  * 0~  r and  (T^  ♦ ...  * c~  p_-L  = m 

one  could  prescribe  a~  arbitrary  analytic  functions  of  p variables 

2 

in  determining  the  conservation  laws,  0"  > 0.  This  means  that  0 for 

example,  could  be  completely  specified  in  a neighborhood  of  U,  which  is 
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2 

clearly  impossible  because  0 satisfies  several  second  order  equations. 
Thus  the  following  important  result  is  clear. 


THEOREM  D;  The  system 


TT(  i , j,k)  = 0 


dTT(  l ,j,k)  = 0 


(3.1*7) 


^nr  * 


i _ 


= 0 


i,j,k  = l,...,p 


is  always  in  involution. 


The  importance  of  this  result  is  that  it  gives  a perfectly  definite 

criterion  for  establishing  the  existence  of  conservation  laws.  One  merely 

needs  to  compute  CT"  and  a~  in  order  to  find  how  arbitrary  the  conservation 
0 -L 

laws  are,  with  no  further  prolongation  necessary. 

When  p = 3 THEOREM  D leads  to  a very  strong  result  since  there  is  at 

most  one  finite  relation,  namely TT « 0.  There  are  then  at  least  eight 
> > 

variables,  n = o,  so  that  n - p “ 5,  unless  some  of  the  variables 

, ^“22’  ^3"*  n°^  aPPear»  an  any  event  n _ P ’ But,  according 

to  the  preceding  theorem,  (T^  = q so  that  q = h.  Hence 

__  > 

(T~  * ” ^ = n - p * 2,  which  implies 


THEOREM  E;  For  p = 3 any  system  of  the  form 


3 

:i.s>)  ^ (* 

<r  =i 


kl  b 

'<r  6X] 


12  3 

~ 9 9 
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satisfying  the  conditions  of  Chapter  I,  § 3,  contains  at  least  a 
two- parameter  family  of  conservation  laws. 


^6.  Several  independent  variables 

The  restrictions  of  Chapter  I,  ^3,  that  the  original  system  should 
contain  only  two  independent  variables,  that  it  should  be  well-determined 
with  distinct  characteristic  roots,  and  that  its  coefficients  should  be 
independent  of  x = { x j , have  been  very  useful,  since  they  permitted  the 
introduction  of  the  normal  form  of  Chapter  I,  § U.  Without  this  normal 
form  the  method  of  Chapter  II  would  have  been  impossible.  However,  it  is 
possible  to  .use  the  methods  of  the  present  chapter  in  a way  which  is 
completely  independent  of  a normal  form  sc  that  the  preceding  restrictions 
become  superfluous.  Here  we  obtain  an  exterior  differential  system  whose 
solutions  correspond  to  conservation  laws  contained  in  the  original  system 
(1.1),  with  no  restrictions  whatsoever  except  that  of  analyticity,  which  is 
required  for  the  Cartan-Kahler  theorem.  No  attempt  will  be  made  to  decide 
when  the  resulting  exterior  differential  system  is  in  involution;  in  order 
to  do  this  it  would  again  be  wise  to  introduce  some  kind  of  normal  form. 
System  (l.l)  is  repeated  here  for  convenience  as 


(3.63) 


-i=1 


ki  ^ 


0, 


k = 1, . . . ,m. 


ki 

Temporarily  we  assume  that  . is  an  analytic  function  of  u 

J 

alone;  this  restriction  will  be  removed  presently. 


-66- 


Suppose  that  there  exist  functions  X ^ such  that  the  equation 


(3.6U) 


n,P,m  v i 

• - 0 


i,  j7k=l 


k i 


has  the  form 


(3.65) 


I 

< 


i=l 


that  is , 

( 3.66) 

Then  clearly 

(3.67) 


n,p  . • 

< ^ - 0. 

. >_i  6xx 


p x i 

du^)  * dd0i  = 0,  i * l,...,n 


by  Poincare's  theorems  so  that  the  necessary  and  sufficient  condition 
for  (3.6U)  to  be  a conservation  law  is  that 


p,m 

(3.68)  d(  ^Xk<^duJ)  =0,  i - l,...,n, 

j,K=l 

/ 

the  sufficiency  following  by  the  converse  of  Poincare's  theorem.  Thus 
the  existence  of  conservation  laws  is  reduced  to  the  problem  of 
prolonging  (3.68)  to  a closed  system  in  involution. 

Now  let  us  consider  an  inhomogeneous  system 


r 


where  o(.  k^  and  oc  k are  analytic  functions  of  x = fx  \ well  a?  u. 
0 


If  there  exist  functions  X , such  that 

K 


(3.70) 


n,P,m  v i 5- 

$ x^fh'  i***  -° 


< k 3 ^ 

i, j,k*l  k=l 


has  the  form  of  a conservation  law  (3.65),  which  in  this  case  becomes 


i,j=i  1=1 


then  clearly 


(3.72)  d ( £ ^ duj  ♦ dx1) 

^ nr 


i = 1, . . . ,n 


that  is , 

P,TO 

(3.73)  d ( ^ Xk-<  t ± dx1  ) - 0,  i = 1, . . . ,n 

J,1<=1 

is  a necessary  and  sufficient  condition  for  the  existence  of  a conservation 
law,  where 


(3.7b) 


i=l  lc=l 


<-v  ^ k 


X v°<- 


i 
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